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Description 

[C The present application claims priority benefit under U.S.C. 119(e) from U.S.' Provisional Application No. 

6C 46,046, filed July 27, 1999, and titled "SYSTEM AND METHOD FOR CONTROL USING QUANTUM SOFT COM- 
PUTING." 

Background of the Invention 
Field of the Invention 

[0002] The present invention relates to programming of quantum computing algorithms, and use of quantum com- 
puting algorithms in connection with control systems based on quantum soft computing. 

Description of the Related Ar t 

[0003] The interplay between mathematics and physics has always been beneficial to both fields of endeavor. The 
calculus was developed by Newton and Leibniz in orderto understand and describe dynamical law of motion of material 
bodies. In general, geometry and physics have had a long and successful symbiotic relationship: classical mechanics 
and Newton's gravity are based on Euclidean geometry, whereas in Einstein's theory of general relativity the basis is 
provided by non-Euclidean, Riemannian geometry (an important insight taken from mathematics into physics). 
Although this link between physics and geometry is still extremely strong, one of the most striking connections today is 
between information theory and quantum physics. There appears to be a trend to make mathematics more "physical." 
[0004] Computation, based on the laws of classical physics, leads to completely different constraints on information 
processing than computation based on quantum mechanics (as first realized by Feymann and Deutsch). Computation 
seems to be the only commodity ever to become exponentially better (e.g., faster) as it gets cheaper. In the past few 
decades, information handling capacity has grown at a rate ten million times faster than that of the human nervous sys- 
tems during the four billion years since life began on Earth. Yet the theory and technology of computing has rested for 
more than 50 years on the Turing-machine model of computation, which leads to many intractable or undecidable prob- 
lems. 

[0005] Quantum computers hold promise for solving such intractable problems, but, unfortunately, there currently 
exist no algorithms for "programming" a quantum computer. Calculation in a quantum computer, like calculation in a 
conventional computer, can be described as a marriage of quantum hardware (the physical embodiment of the comput- 
ing machine itself, such as quantum gates and the like), and quantum software (the computing algorithm implemented 
by the hardware to perform the calculation). To date, quantum software algorithms, such as Shor's algorithm, used to 
solve problems on a quantum computer have been developed on an ad hoc basis without any real structure or program- 
ming methodology. 

[0006] This situation is somewhat analogous to attempting to design a conventional logic circuit without the use of 
a Karnough map. A logic designer, given a set of inputs and corresponding desired outputs, could design a complicated 
logic circuit using NAND gates without the use of a Karnough map. However, the unfortunate designer would be forced 
to design the logic circuit more or less by intuition, trial, and error. The Karnough map provides a structure and an algo- 
rithm for manipulating logical operations (AND, OR, etc.) in a manner that allows a designer to quickly design a logic 
circuit that will perform a desired logic calculation. 

[0007] The lack of a programming or design methodology for quantum computers severely limits the usefulness of 
the quantum computer. Moreover, it limits the usefulness of the quantum principles, such as superposition, entangle- 
ment and interference, that give rise to the quantum logic used in quantum computations. These quantum principles 
suggest, or lend themselves, to problem-solving methods that are not used in conventional computers. 
[0008] These quantum principles can be used with conventional computers in much the same way that genetic prin- 
ciples of evolution are used in genetic optimizers today. Nature, through the process of evolution, has devised a useful 
method for optimizing large-scale nonlinear systems. A genetic optimizer running on a computer efficiently solves many 
previously difficult optimization problems by simulating the process of natural evolution. Nature also uses the principles 
of quantum mechanics to solve problems, including optimization-type problems, searching-type problems, selection- 
type problems, etc. through the use of quantum logic. However, the quantum principles, and quantum logic, have not 
been used with conventional computers because no method existed for programming an algorithm using the quantum 
logic. 
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Summary of the Invention 

[0009] The present invention solves these and other problems by providing a methodology and an algorithm for pro- 
gramming an algorithm to solve a problem using quantum logic. The quantum logic program can be "run" on a quantum 
computer. The algorithm can also be "run" on a non-quantum computer by using the non-quantum computerto simulate 
a quantum computer. This allows the concepts, features, and principles of quantum computing, such as superposition, 
entanglement, quantum interference, and the like (and the massive parallelism enabled by these principles) to be used 
to advantage in non-quantum computers without the need to develop quantum computer hardware. 
[0010] In one embodiment, the quantum programming method is used to with a genetic search algorithm in a con- 
trol system. A conventional genetic search algorithm searches for an optimal solution in a single space. The quantum 
search algorithm provides global searching for an optimum solution among many spaces. 

[0011] In one embodiment, an algorithm design for quantum soft computing is designed by encoding an input func- 
tion f into a unitary matrix operator U F . The operator U F is embedded into a quantum gate G, where G is a unitary 
'matrix. The gate G is applied to an initial canonical basis vector to produce a basis vector. The basis vector is meas- 
ured. These steps are repeated several times as necessary to generate a set of measured basis vectors. The measured 
basis vectors are decoded and translated into' an output vector. 

[001 2] In one embodiment, the encoding into U F includes transforming a map table of f into an injective function F, 
transforming the map table of F into a map table for.U F and transforming the map table for U F into U F . 
[0013] In one embodiment, the Shannon entropy of the basis vectors is minimized. 

[0014] In one embodiment, an intelligent control system having a quantum search algorithm to reduce Shannon 
entropy includes a genetic optimizer to construct local solutions using a fitness function configured to minimize a rate 
of entropy production of a controlled plant. A quantum search algorithm is used to search the local solutions to find a 
global solution using a fitness function configured to minimize Shannon entropy. 

[0015] In one embodiment, global solution includes weights for a fuzzy neural network. In one embodiment, the 
fuzzy neural network is configured to train a fuzzy controller, and the fuzzy controller provides control weights to a pro- 
portional-integral-differential (PID) controller. The PID) controller controls a plant. 

[0016] In one embodiment, a quantum search algorithm is evolved according to a fitness function selected to min- 
imize Shannon entropy. 

[0017] In one embodiment, quantum search algorithm is evolved by minimizing Heisenberg uncertainty and mini- 
mizing Shannon entropy. 

[0018] In one embodiment, a quantum search evolves by applying an entanglement operator to create a plurality of 
correlated state vectors from a plurality of input state vectors and applying an interference operator to the correlated 
state vectors to generate an intelligent state vector, where the intelligent state vector has less classical entropy than the 
correlated state vectors. 

[0019] In one embodiment, global optimization to improve a quality of a sub-optimal solution is accomplished by 
applying a first transformation to an initial state to produce a coherent superposition of basis states. A second transfor- 
mation is applied, to the coherent superposition using a reversible transformation to produce coherent output state. A 
third transformation is applied to the coherent output states to produce an interference of output states and a global 
solution is selected from the interference of output states. In one embodiment, the first transformation is a Hadamard 
rotation. In one embodiment, each of the basis states is represented using qubits. In one embodiment, the second 
transformation is a solution to Shrodinger's equation. In one embodiment, the third transformation is a quantum fast 
Fourier transform. In one embodiment, selecting is made to find a maximum probability. In one embodiment, the super- 
position of input states includes a collection of local solutions to a global fitness function. 

Brief Description of the Drawings 

[0020] 

Figure 1A is a block diagram showing a control' system that uses soft computing in combination with a genetic 
search algorithm. 

Figure 1 B is a block diagram showing the process of the genetic search. 

Figure 2 is a block diagram showing a control system that uses soft computing in combination with a quantum 
search algorithm. 

Figure 3 is a block diagram showing the structure of a classical genetic algorithm and a quantum searching algo- 
rithms for global optimization. 

Figure 4 is a block diagram showing a general structure of a quantum search algorithm. 
Figure 5 is a block diagram showing a quantum network for a quantum search algorithm. 
Figure 6 is a block diagram of the quantum searching algorithm. 
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Figure 7 illustrates a gate approach for simulation of quantum algorithms using classical computers. 

Figure 8 illustrates a programming diagram for quantum algorithms. 

Figure 9 illustrates the structure of the quantum encoder shown in Figure 8. 

Figure 10 illustrates the structure of the quantum block shown in Figure 8. 
5 Figure 1 1 illustrates an example of a quantum circuit. 

Figure 12A illustrates a quantum circuit for a tensor product transformation. 

Figure 12B illustrates a quantum circuit for a dot product transformation. 

Figure 12C illustrates a quantum circuit for an identity transformation. 

Figure 12D illustrates a quantum circuit for propagation. 
70 Figure 1 2E illustrates a quantum circuit for iteration. 

Figure 12F illustrates a quantum circuit for input/output. 

Figure 13 illustrates a first representation of a quantum circuit for the Deutsch-Jozsa's Quantum Gate. 

Figure 14 illustrates a second representation of a quantum circuit for the Deutsch-Jozsa's Quantum Gate. 

Figure 15 illustrates a circuit for the Constant Function with Value 1 - first circuit. 
15 Figure 1 6 illustrates a circuit for the Constant Function with Value 1 - second circuit. 

Figure 17 illustrates a circuit for the Constant Function with Value 1 - third circuit. 

Figure 18 illustrates a circuit for the Constant Function with Value 1 - fourth circuit. 

Figure 19 illustrates a circuit for the Constant Function with Value 1 - fifth circuit. 

Figure 20 illustrates a constant Function with value 0. 
20 Figure 21 illustrates the evolution of the DJ quantum gate. 

Figure 22 illustrates the final DJ quantum gate. 

Figure 23 illustrates the quantum circuit of Graver's Quantum Gate. 

Figure 24 illustrates the final circuit for Grover's Quantum Gate. 

Figure 25 is a Wenn diagram that illustrates entropy and mutual information in a quantum system. 
25 Figure 26 is a general circuit diagram for a quantum block. 

Figure 27 illustrates information analysis of the Deutsch-Jozsa algorithm for a first operator U F . 

Figure 28 illustrates information analysis of the Deutsch-Jozsa algorithm for a second operator U F . 

Figure 29 illustrates information analysis of the Deutsch-Jozsa algorithm for a third operator U F . 

Figure 30 illustrates information analysis of the Deutsch-Jozsa algorithm for linear functions. 
30 Figure 31 illustrates information analysis of the Deutsch-Jozsa algorithm for non-linear balanced functions. 

Figure 32 illustrates the quantum information structure of the Deutsch-Jozsa input space. 

Figure 33 shows the quantum information structure of the Shor input space. 

Figure 34 illustrates information analysis of the Shor algorithm for a first operator U F . 

Figure 35 illustrates information analysis of the Shor algorithm for a second operator U F . 
35 Figure 36 illustrates information analysis of the Grover algorithm for a general iteration. 

Figure 37 illustrates information analysis of the Grover algorithm for a first iteration. 

Figure 38 illustrates information analysis of the Grover algorithm for a second iteration. 



40 

Introduction 

[0021] Classical control theory is based on the assumption that all controlled plants can be approximated as linear 
systems. Unfortunately, this assumption is rarely true in the real world. Most plants are highly non-linear, and often have 
45 no simple control algorithm. To meet this need for nonlinear control, current control systems have been developed that 
use soft computing concepts such as genetic analyzer, fuzzy neural networks and the like. The control system evolves 
(changes) over time to adapt to changes in the plant, the operating environment, or both. 

[0022] Figure 1 A shows a control system 1 00 based on soft computing to control a plant 1 04. An entropy production 
calculator 1 06 calculates an entropy production rate for the plant 1 04. An output from the entropy production calculator 
so 1 06 is provided to an input of a Genetic Analyzer (GA) 1 07 (the GA can also be called a genetic optimizer). The GA 1 07 
searches for a set of control weights that produce minimum entropy production. The weights are provided to a fuzzy 
neural network (FNN) 1 08. An output of the FNN is provided to a fuzzy controller 1 09. An output of the fuzzy controller 
1 09 is a set of gain schedules for a conventional proportional-integral-differential (PID) controller 1 03 that controls the 
plant 1 04. 

55 [0023] For soft computing systems' based on a genetic analyzer (GA), there is very often no real control law in the 
classic control sense out rather, control can be based on a physical control law such as minimum entropy production. 
This allows robust control because the GA 1 07, combined with feedback, guarantees robustness. However, robust con- 
trol is not necessarily optimal control. 
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[0024] The GA 107 attempts to find a global optimum solution for a given solution space. Any random disturbance 
of the plant can "kick" the GA 107 into a different solution space. Thus for example, if the controlled plant is a suspen- 
sion system of a car, a change in road conditions can force the GA into a different solution space. 
[0025] Genetic algorithms are global search algorithms based on the mechanics of natural genetics and natural 

5 selection. In the genetic search, each design variable is represented by a finite length binary string and these finite 
binary strings are connected in a head-to-tail mannerto form a single binary string. Possible solutions are coded to rep- 
resented by a population of binary strings. Genetic transformations analogous to biological reproduction and evolution 
are subsequently used to improve and vary the coded solutions. Usually, three principal operators i.e., reproduction, 
crossover, and mutation are used in the genetic search. 

10 [0026] The reproduction process is one that biases the search toward producing more fit members in the population 
and eliminating the less fit ones. Hence, a fitness value is first assigned to each string in the population. One simple 
approach to select members from a initial population to participate in the reproduction is to assign each member a prob- 
ability of selection on the basis of its fitness value. A new population pool of the same size as the original is then created 
with a higher average fitness value. 

15 [0027] The process of reproduction simply results in more copies of the dominant design to be present in the pop- 
ulation. The crossover process allows for an exchange of design characteristics among members of the population pool 
with the intent of improving the fitness of the next generation. Crossover is executed by selecting strings of two mating 
parents, randomly choosing two sites on the strings, and swapping strings of O's and 1 's between these chosen sites. 
[0028] Mutation safeguards the genetic search process from a premature loss of valuable genetic material during 

20 reproduction and crossover. The process of mutation is simply to choose a few members from the population pool 
according to the probability of mutation and to switch a 0 to 1 or vice versa at a randomly selected mutation site on the 
selected string. 

[0029] A simple schematic illustration of the process of genetic search (as used by the GA 1 07) is shown in Figure 
1B. A population of strings is first transformed into decimal codes and then sent into the fitness process 160 forcom- 

25 puting the fitness of all strings. A biased roulette wheel 1 06 where each string has roulette wheel slot sized in proportion 
to its fitness is created. A spinning of the weighted roulette wheel yields the reproduction candidate. In this way, a higher 
fitness of strings has a higher number of offspring in the succeeding generation. Once a string has been selected for 
reproduction, a replica of the string based on its fitness is created and then entered into a mating pool 1 62 for waiting 
the further genetic operations. After reproduction, a new population of stings is generated trough the evolution proc- 

30 esses of crossover 163 and mutation 164. Finally, the whole genetic process, as mentioned above, is repeated again 
and again until an optimal solution is found. 

[0030] In the foregoing discussion, the mechanics of the genetic search are simple, however, there are some key 
^differences from traditional methods which contribute to the strengths of the approach. Genetic algorithms work on 
function evaluations alone and do not require function derivatives. While derivatives contribute to a faster convergence 

35 towards the optimum, they may also direct the search towards a local optimum. Furthermore, since the search pro- 
ceeds from several points in the design space to another such set of design points, the method has a better probability 
of locating a global minimum as opposed to those schemes that proceed from one point to another. In addition, genetic 
algorithms work on a coding design variables rather than the variables themselves. This allows for an extension of these 
algorithms to design space consisting of a mix of continuous, discrete, and integer variables. 

40 [0031] As stated above the GA 1 07 searches for a global optimum in a single solution space. It is desirable, how- 
ever, to search for a global optimum in multiple solution spaces to find a "universal" optimum. A quantum search algo- 
rithm provides the ability to search multiple spaces simultaneously (as described below). Figure 2 shows a modified 
version of the control system wherein a Quantum Search Algorithm (QSA) 202 is interposed between the GA 107 and 
the FNN 108. The QSA 108 searches several solution spaces, simultaneously, to find a universal optimum, that is, a 

45 solution that is optimal considering all solution spaces. 

[0032] Quantum Algorithms are global random searching algorithms based on the quantum mechanics principles, 
laws, and quantum effects. In the quantum search, each design variable is represented by a finite linear superposition 
of classical initial states, with a sequence of elementary unitary steps manipulate the initial quantum state li> (for the 
input) such that a measurement of the final state of the system yields the correct output. It begins with elementary clas- 

50 sical preprocessing, and then it applies the following quantum experiment: starting in a initial superposition of all possi- 
ble states, it computes a classical function, applies a quantum Fourier transform, and finally, performs a measurement. 
Depending on the outcome, it may carry out one or more similar quantum experiments, or complete the computation 
with some classical post-processing. Usually, three principal operators are used in the quantum searching algorithm. 
These operators are linear superposition (coherent states), entanglement, and interference. 

55 [0033] Linear superposition is closely related to the familiar mathematical principle of linear combinations of vec- 
tors. Quantum systems are described by a wave function y that exists in a Hilbert space. The Hilbert space has a set 
of states l<jty, that form a basis, and the system is described by a quantum state: 
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!\|>> is said to be in a linear superposition of the basis states i<^, and in the general case, the coefficients c, may be 
complex. Use made here of the Dirac bracket notation, where the ke I > is analogous to a column vector, and the 
bra ( I is analogous to the complex conjugate transpose of the ket. 

[0034] In quantum mechanics, the Hilbert space and its basis have a physical interpretation, and this leads directly 
to perhaps the most counterintuitive aspect of the theory. The counter intuition is that (at the microscopic level), the 
state of the system is described by the wave function, that is, as a linear superposition of all basis states (i.e., in some 
sense, the system is in all basis states at once). However, at the macroscopic or classical level, the system can be in 
only a single basis state. For example, at the quantum level, an electron can be in a superposition of many different 
energies; however, in the classical realm, this cannot be. This is similar to saying that during a coin toss, while the coin 
is in the air, the coin is in a superposition of both of its basis states (heads and tails). When the coin lands, the wave 
function "collapses" into one of the classical states, either heads ortails. While the coin is in the air, one cannot say how 
it will land, one can only assign probabilities. The quantum mechanical wave function x|/ does the same thing. The wave 
function \|/ assigns probabilities that the system will "land" in each of its basis states. 

[0035] Coherence and decoherence are closely related to the idea of a linear superposition. A quantum system is 
said to be coherent if it is in a linear superposition of its basis states. A result of quantum mechanics is that if a system 
is in a linear superposition of states interacts with its environment, the superposition is destroyed. This loss of coher- 
ence is called decoherence and is governed by the wave function. The coefficients q are called probability amplitudes 
and lc,l 2 gives the probability of l\|/> collapsing into the state ty> upon the occurrence of an interaction with the environ- 
ment that causes decoherence (i.e. a measurement). The wave function describes a real physical system that must col- 
lapse into exactly one basis state. Therefore, the probabilities governed by the amplitudes c ; must sum to unity. This 
necessary constraint is expressed as the unitary condition 

£lc,.| 2 =1 

[0036] In the Dirac notation, the probability that a quantum state will collapse into an eigenstate (basis state) is writ- 
ten l^/l^l 2 and is analogous to the dot product (projection) of two vectors. Consider, for example, a discrete physical 
variable, called spin. The simplest spin systems is a two-state system called a spin- 1/2 system, whose basis states are 
represented as spin-up and spin-down. In this simple system the wave function is a distribution over two values (up and 
down) and a coherent state is a linear superposition of spin up and spin down. As long as the system maintains its 
quantum coherence, it cannot be said to be either spin up or spin down. It is in some sense both at once. Classically, 
of course, it must be one or the other, and when this system undergoes decoherence, the result is either sin up or spin 
down. 

[0037] Operators on a Hilbert space describe how one wave function is changed into another. They will typically be 
denoted as matrices acting on vectors. Using operators, an eigenvalue equation can be written 



where a,- are the eigenvalues of the operator A. The solutions to such an equation are called eigenstates and can 
be used to construct the basis of a Hilbert space. In the quantum formalism, all properties are represented as operators 
whose eigenstates are the basis for the Hilbert space associated with that property and whose eigenvalues are the 
quantum allowed values for the at property. It is important to note that operators in quantum mechanics are linear oper- 
ators and further that they are unitary so that 

A*A = AA*=I, 

where I is the identity operator, and A* is the complex conjugate transpose, (also known as the adjoint) of A. 
[0038] Interference is a familiar wave phenomenon. Wave peaks that are in phase will constructively add, while 
those that are out of phase interfere destructively. This is a phenomenon common to all kinds of wave mechanics from 
water waves to light. The well known double slit experiment demonstrates empirically that at the quantum level, inter- 
ference also applies to the probability waves of quantum mechanics. 
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[0039] Entanglement is the potential for quantum states to exhibit correlations that cannot be accounted for classi- 
cally. Form a computational standpoint, entanglement seems intuitive enough, it is simply the fact that correlations can 
exist between different quantum systems (e.g. different particles). For example, if one quantum particle is the in spin- 
up state, then another will be in the spin-down state. What makes it. so powerful is the fact that since quantum states 
5 exist as superpositions, these correlations exist in the superposition as well. When the superposition is destroyed, the 
proper correlation is instantly communicated between the particles, and it is this communication that is the crux of 
entanglement. The communication occurs instantaneously, even if the particles are separated by a large distance. Ein- 
stein called this "spooky action at a distance." 

[0040] From a computational standpoint, quantum states that are superpositions of only basis states that are max- 
10 imally far apart in terms of Hamming distance are those states with the greatest entanglement. Moreover, while inter- 
ference is a quantum property that has a classical cousin; entanglement is a completely quantum phenomenon for 
which there is no classical analog. 

[0041] Quantum networks are one of the several models of quantum computation. Others include quantum Turing 
machines, and quantum cellular automa. In the quantum networks model, each unitary operator is modeled as a quan- 
15 turn logic gate that affects one, two, or more quantum bits (qubits). Qubits exist in a superposition of states, thus, quan- 
tum logic gates operate on qubits by operating on all states in the superposition simultaneously. This results in quantum 
parallelism. (The term quantum logic gate is simply a schematic way to represent the time evolution of a quantum sys- 
tem. The term "gates" is not meant to imply that quantum computation can be physically realized in a manner similar to 
classical logic networks.) 

20 [0042] Because of the entanglement (quantum correlations) between the quantum particles, the state of the system 
cannot be specified simply by describing the state of each particle. Instead, the state of the particles is a complicated 
superposition of all 2" basis states, so 2" complex coefficients are needed in order to describe it. The exponentially 
large size of a Hilbert space is an ingredient in quantum computations. It is tempting to think that exponential parallelism 
implies exponential computational power, but this is not the case. In fact, classical computations can be viewed as hav- 

25 ing exponential parallelism as well, the problem lies in the question of how to extract the exponential information out of 
the system. In quantum computation, in order to extract quantum information, one has to observe the system. The 
measurement process causes the famous collapse of the wave function. In a nutshell, the means that after the meas- 
urement, the state is projected onto one of the exponentially many possible states, so that the exponential amount of 
information which has been computed is lost. 

so [0043] In orderto take advantage of exponential parallelism, one needs to combine it with interference. Interference 
allows the exponentially many computations done in parallel to cancel each other, just like destructive interference of 
waves or light. The goal is to arrange the cancellation such that only the computations which are of interest remain, and 
: all the rest cancel out. The combination of exponential parallelism and interference is what makes quantum computa- 
tions powerful and plays an important role in quantum algorithms. The Fourier transform indeed manifests interference 

35 and exponentiality. 

[0044] In classical computations and in digital electronics, one deals with sequences of elementary operations 
(operations such as ANDS, OR, and NOT). These sequences are used to manipulate an array of classical bits. The 
operations are elementary in the sense that they act on only a few bits (one or two) at a time. It is sometimes convenient 
to refer to sequences as products and the operations as operators, matrices, instructions, steps, or gates. In quantum 
40 computation, one also deals with sequences of elementary operations (such as controlled-NOT and qubit rotations), but 
for manipulating qubits instead of classical bits. Quantum sequences of elementary operations are often represented 
graphically by qubit circuits. In quantum computations, one often knows the unitary operator U that describes the evo- 
lution of an array of qubits. 

[0045] The exponential parallelism inherent in quantum systems allows quantum searching algorithms to provide 
45 for the support of decision-making process in control systems. Figure 3 compares the structures of the GA and QSA 
algorithms. As shown in Figure 3, in the GA search, a solution space 301 leads to an initial position (input) 302. The 
initial position 302 is coded into binary strings using a binary coding scheme 31 0. GA operators such as selection 303, 
crossover 304, and mutation 305 are applied to the coded strings to generate a population. Through a fitness function 
306 (such as a fitness function based on minimum entropy production or some other desirable property) a global opti- 
50 mum for the space 301 is found. 

[0046] By contrast, in the QSA shown in Figure 3, a group of N solution spaces 350 are used to create an initial 
position (input) 351 . Quantum operators such as superposition 352, entanglement 353, and interference 354 operate 
on the initial position to produce a measurement. Superposition is created using a Hadamard transformation 361 (a 
one-bit operation). Entanglement is created through a Controlled-NOT operation 362 (a two-bit operation). Interference 
55 is created through a Quantum Fourier Transform (QFT) 363. Using the quantum operators, a universal optimum for cov- 
ering all the spaces in the group 350 is found. 

[0047] Thus, the classical process of selection is loosely analogous to the quantum process of creating a superpo- 
sition. The classical process of crossover is loosely analogous to the quantum process of entanglement. The classical 
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process of mutation is loosely analogous to the quantum process of interference. 

[0048] Figure 4 shows a general structure of a QSA (such as the QSA 202) having a conceptual level 400, a struc- 
ture level 401 , a hardware level 402, and a software level 403. At the conceptual level 400, an initial state 41 0 is pro- 
vided to a process block 420 that creates a superposition of states. The superposition of states is provided to a process 
s block 430 that provides a unitaary operator U f to the superposition. An output of the process block 430 is provided to a 
solution block 440 that computes an interference of solutions. An output of the solution block 440 is provided to a obser- 
vation/measurement block 460. 

[0049] At the structure level, an input is coded as a series of quantum bits (qubits) that are prepared in an initial 
state (e.g., a logical zero state) and provided to a Hadamard Transformation Matrix 421 to produce the superposition. 
10 The superposition from the matrix 421 is provided to the operator U f where U f is a solution to the Schrodinger Equation 
in a process block 431. An output .from the process block, being a solution to the Schrodinger equation, is provided to 
a Quantum Fast Fourier Transform (QFFT) 441 to provide the interference. The output of the QFFT 441 is provided to 
a Transformation Matrix 451 . An output of the Transformation matrix 451 is provided as a solution of maximum proba- 
bility amplitude 461 . 

15 [0050] At the hardware level, the superposition 420 is produced by rotation gates 422, the operator U f is imple- 
mented as a sequence of elementary gate operations 432, the QFFT 441 is implemented as a sequence of Hadmarad 
ana F-srmutation (P) operator gates, and the Transformation Matrix 451 is implemented using rotation gates 452. 
[0051] Figure 5 illustrates the QSA as an architecture involving the sequence from an initial state, through the cre- 
ation of a superposition. Entanglement is applied to the superposition using quantum parallelism inherent in a coherent 

20 quantum system with entangled states. The parallelism collapses when interference is introduced to produce a super- 
position of solutions through the QFFT. Figure 5 illustrates these processes by comparing the classical double slit 
experiment to logical quantum operations and to the quantum search operations. 

[0052] In the classical double slit, a source 501 creates a particle having an initial superposition of states. This is 
analogous to the quantum algorithm operation of applying a Hadamard (rotation gates) to a qubit initial to an eigenstate. 
25 [0053] Returning to the double slit, entanglement is produced with the particle passes through slits 502. This cor- 
responds to the process of operating on the superposition using the unitary operator U { . 

[0054] Again returning to the double slit, interference is produced when the entangled particles strike a photo- 
graphic film placed behind the slits to produce an interference pattern ( a superposition of solutions). This corresponds 
to the QFFT. 

30 [0055] Finally, selection of the desired solution corresponds to choosing the larges probability from the QFFT (that 
is, the brightest line produced on the film). 

[0056] Figure 6 shows use of the QSA in connection with a GA 605. In Figure 6, a generator of initial states 604 
works in connection with the GA 605 and, optionally a fuzzy neural network 603, to produce a set of initial states. The 
initial states are provided to a Hadamard transform 602 to produce a superposition of classical states 601 . The super- 
35 position of classical states is provided to a process block 606 that introduces entanglement through the use of operators 
such as a controlled NOT. An output from the process block 606 is provided to an interference block 607 that computes 
an interference of the entangled states using a QFFT. An output from the interference block 607 is provided to a meas- 
urement/observation block 608 which selects a desired solution from the superposition of solutions computed by the 
block 607. 

40 [0057] An output from the observation/measurement block 608 is provided to a decision block 609. The decision 
block 609 makes decisions such as inputs for the generator of initial states 604 and, optionally, a new fitness function 
for the GA 605. The decision block 609 can also provide data to, and receive data from a decoding block. 61 0. The 
decoding block 610 can communicate with sensors, other control systems, users, etc. 

[0058] The basis of quantum computation is obtained from the laws of quantum information theory wherein infor- 
ms mation is something that encoded in the state of a physical system, and a computation is something that can be carried 
out on an actual physical realizable device. 

[0059] Two examples can help demonstrate the effects of quantum constrained approach to develop a correctness 
model of physical objects: First, comparing a classical Shannon data transmission channel C c , to quantum data trans- 
mission channel C q shows that: 

50 
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where P is an input signal of a channel with the power P and N is a power of a noise. For N -> 0, lim C c -> 
lim C q -» In P and we have the correctness model of information transmission channel with noise. 
[0060] Second, consider identification of a mathematical model of a control object with a quantum constraint. For a 
classical object, [x,y] = 0 (commutative case), 

y(t)=jx(x)k(t-%)dx (2.1) 



[0061] Whereas, for an object based on quantum correction of model identification, 

[x,y] = ih 

(non-commutative case) 

At) = ^x(t)+)x(T)k(t-T)dT. ' (2.2) 



In the quantum case, non-correctness of a Fredholm integral equation of 1st order maps into correctness of a Fredholm 
integral equation of 2nd order. This means that it is possible with the small error in measurements of an output signal 
25 y(t) in (2.2) to do the correctness identification of the input signal x(t) or a structure of dynamical system k(t) also with 
small error. In the classical case (2.1 ), a non-correctness model identification and a small error in measurement output 
signal gives a large error in the identification signal. 

[0062] Quantum computing is based on the principles of quantum superposition, quantum interference and quan- 
tum entanglement. The quantum state at all times has components corresponding to some or all of the possible clas- 

30 sical states. This quantum effect is known as a superposition state. A computer built upon quantum rules can process 
different inputs using massive parallelism and produce a superposition of outputs. In this case a quantum computer is 
a physical machine that can accept input states which represent a coherent superposition of many different possible 
inputs and subsequently evolve them into a corresponding superposition of outputs. Quantum entanglement allows one 
to encode data into non-trivial multi-particle superpositions of some pre-selected basis states, and quantum 

35 interference (which is a dynamic process) allows one to evolve initial quantum states (as inputs) into final states (as out- 
puts) modifying intermediate multi-particle superpositions in some prescribed way. Quantum computers use the quan- 
tum interference of different computational paths to enhance correct outcomes and suppress erroneous outcomes of 
computations. A common pattern underpinning quantum algorithms can be identified when quantum computation is 
viewed as multi-particle interference. The multi-particle interference (unlike single-particle interference) does not any 

40 classical analogue and can be viewed as inherently quantum process. 

[0063] Classical computers, although based on quantum physics, are not fully quantum, they do not use " quantum- 
ness* of matter at the information-theoretic level, where it really matters. Namely, in a classical computer information is 
recorded in macroscopic, two level systems. Wires conducting electrical current in computers can be in two basic 
states: when there is no current flowing through, representing a logical "0", or else when there is some current flowing 

45 through, representing a logical "1". These two states form a bit of information. All computation is based on logical 
manipulation of bits through logical gates acting on wires representing these bits. However, as described herein, a 
quantum computer can be simulated on a classical computer. 

[0064] Instead of wires and currents, a quantum computer uses, for example, two electronic states of an atom to 
record information. The two quantum states can be, for example, the ground state, I0>, and the excited state, I1> (using 
so Dirac notation). Since an atom obeys the laws of quantum mechanics, the most general electronic state is a superpo- 
sition of the two basic states given by hjf = 0,10) + c 2 M>, and called the quantum-bit or qubit, for short. In addition 
to 0 and 1 states, a qubit has, so to speak, all the states "in between". With two classical bits, there are four possibilities: 
00, 01, 10, 11. This should be contrasted with two qubits which are in general in a state of the form: 
IVi>= c-,IOO> + c 2 I01>+ c 3 110)+ c 4 IH>. If for example c., = c 4 = 0 and 
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n one has the famous Einstein-Podolski-Rosen (EPR) state: 



[0065] Two qubits in this state display a degree of correlation impossible in classical physics and hence violate the 
Bell inequality, which is satisfied by all local (i.e. classical) states. This phenomenon is called entanglement and is at 
the root of the success of quantum computing. The exploitation of a number of entangled qubits can lead to a consid- 
10 erable computational speed-up in a quantum computer over its classical counterpart. Thus, what distinguishes classical 
and quantum computing is how the information is encoded and manipulated, i.e. what plays a crucial role is whether 
the logical basis is the classical (Boolean logic) or the quantum logic. 

[0066] A quantum computer can solve problems such as factoring and database search very efficiently. Quantum 
searching algorithms (QSA) can be used to solve problems that are not resolvable using classical algorithms. In one 
15 embodiment, quantum mechanical algorithms are combined efficient database algorithms that make jse of specific 
properties of the database. 



General Structure of Quantum Algorithms 



20 [0067] In one embodiment, a Quantum Algorithm, written as a Quantum Circuit, is translated into the corresponding 
Programmable Quantum Gate. This gate is represented as a matrix operator such that, when it is applied to the vector 
representation of the quantum register input state, the produced result is the vector representation of the a register out- 
put state, as shown in Figure 7 

[0068] Quantum computation is based on three operators that act on quantum coherent states. The operators are: 
25 superposition, entanglement and interference. The coherent states are described as those solutions of the correspond- 
ing Schrodinger equation that represent the evolution states with minimum uncertainty (according to Heisenberg, they 
are those quantum states with "maximum classical properties"). The Hadamard Transform creates the superposition 
from classical states, and quantum operators as CNOT create robust entangled states. The Quantum Fast Fourier 
Transform provides interference. The efficient implementation of a number of operations for quantum computation 
30 include controlled phase adjustment of the amplitudes in superposition, permutation, approximation of the transforma- 
tions, and generalizations of the phase adjustments to block-matrix transformations. These operations generalize those 
used in quantum search algorithms that are realized on a classical computer. The following sections demonstrate the 
application of this approach to the general case of simulation on classical computers, followed by examples based on 
the Deutsch-Jozsa algorithm and the Graver algorithm. 
35 [0069] The problems solved by a quantum algorithms can be so stated as: 



given a function f:{0,1}"->{0,1} m find a certain property of f. 



[0070] The structure of a quantum algorithm is outlined, with a high level representation, in Figure 8. In Figure 8, 
40 an input 801 , describing the function f, is provided to an encoder 802. An output of the encoder 802 is an operator U f 
803. The operator 803 is provided to a quantum block 804. An output of the quantum block 804 is a set of basis vectors 
805. The basis vectors 805 are provided to a decoder 806. An answer 807 is provided as an output of the decoder 806. 
The input 801 and output 807 make up a binary string layer. The encoder 802 and decoder 806 make up a map table 
and interpretation space level. The operator 803, the quantum block 804, and the basis vectors 805 exist in a complex 
45 Hilbert space. 

[0071] The input of a quantum algorithm is a function f that maps binary strings into binary strings. This function is 
represented as a map table, defining a mapping for every string its image. The function f is firstly encoded into a unitary 
matrix operator U F that depends on the properties of /. In some sense, the unitary operator U F calculates f when its 
input and output strings are encoded into canonical basis vectors of a complex Hilbert Space. U F maps the vector code 
so of every string into the vector code of its image by /. A squared matrix U F on the complex field is unitary if and only if 
its inverse matrix coincides with its conjugate transpose: U F ~ 1 =l/ F . A unitary matrix is always reversible and pre- 
serves the norm of vectors. 

[0072] Once the matrix operator U F has been generated, it is embedded into a quantum gate G, a unitary matrix 
whose structure depends on the form of matrix U F and on the problem to be solved. The quantum gate is the heart of 
55 a quantum algorithm. The quantum gate acts on an initial canonical basis vector to generate a complex linear combi- 
nation (a superposition) of basis vectors as an output. This superposition contains all of the information needed to 
answer the initial problem. 

[0073] After this superposition has been created, measurement takes place in order to extract the information. In 
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quantum mechanics, measurement is a non-deterministic operation that produces asoutput only one of the basis vec- 
tors in the entering superposition. The probability of every basis vector of being the output of measurement depends on 
its complex coefficient (probability amplitude) in the entering complex linear combination. 

[0074] The segmental action of the quantum gate and of measurement is provided by the quantum block 804. The 
s quantum block 804 is repeated k times in orderto produce a collection of k basis vectors. Since measurement is a non- 
deterministic operation, these basic vectors will not necessarily be identical, and each one will encode a piece of the 
information needed to solve the problem. 

[0075] The last part of the algorithm includes the interpretation of the collected basis vectors in order to get the 
answer to the initial problem with a certain probability. 
w [0076] The behaviour of the encoder block 802 is illustrated in Figure 9. As shown in Figure 9, the function f is 
encoded into matrix U F in three steps. First, the map table of function f:{0,1}"->{0,1} m is transformed into the map table 
of the injective function F:{0,1} n+m ->{0,1} n+m such that: 

F (* 0 . - x„.i, y 0 . ••> y m -i)= (* 0 . - x n-v f (x 0 , ... *n-i)®(y 0 . ... y m -i)) 

[0077] The injective function comes from the requirement that U F is unitary. A unitary operator is reversible, so it 
cannot map two different inputs in the same output. Since U F is the matrix representation of F, F is supposed to be 
injective. To directly use the matrix representation of function f, would give a non-unitary matrix, since f could be non- 
injective. So, injectivity is fulfilled by increasing the number of bits and considering the function F instead of the function 
20 f. The function f can always be calculated from F by putting (y 0> .., y m . 1 )=(0,..,0) in the input string and reading the 
last m values of the output string. 

[0078] In the second step in Figure 9, the function F map table is transformed into a U F map table, according to the 
following constraint: 

25 vse{0,1} n+m : U F [x(s)]=z[F(s)} 

[0079] The code map % : {0,1} n+m -> C 2n+m (C 2n * m is the target Complex Hilbert Space) is such that: 

«■(*, .-. *„_ ) - K»o )® •• « 4w. ) - 1 *«• . ) 

35 

[0080] The code x maps bit values into complex vectors of dimension two belonging to the canonical basis of C 2 . 
Using the tensor product, x maps the general state of a binary string of dimension n into a vector of dimension 2", 
reducing this state to the joint state of the n bits composing the register. Every bit state is transformed into the corre- 
40 sponding 2-dimesional basis vector, and then the string state is mapped into the corresponding 2"-dimesional basis 
vector by composing all bit-vectors through the tensor product. In this sense tensor product is the vector counterpart of 
state conjunction. 

[0081 ] Finally, in step three of Figure 9, the U F map table is transformed into U F using the following transformation 
rule: 

45 

lUJ^I** U F \j) = \,) 

[0082] This rule can be understood by considering vectors l/> and \j> as column vectors. Associating these vectors 
to the canonical basis, U F defines a permutation map of the identity matrix rows. In general, row ly> is mapped into row 
so \i>. 

[0083] Figure 1 0 illustrates the operation of the quantum block 804. The heart of the quantum block 804 is the quan- 
tum gate 1 002, which depends on the properties of matrix U F . 

[0084] The matrix operator U F is the output of the encoder block 802, and it is the input for the quantum block 804. 
In the quantum block 804, the matrix operator U F is first embedded into a more complex gate: the quantum gate G. The 
55 unitary matrix G is applied k times to an initial canonical basis vector l/> of dimension 2 n+m . Every time, the resulting 
complex superposition (3I0..01 ..1 > of basis vectors is measured, producing one basis vector lx,> as a result. The meas- 
ured basis vectors {\x,>,..,\x k >} are collected together. This collection is the output of the quantum block 804. The "intel- 
ligence" of this algorithm lies in the ability to build a quantum gate that is able to extract the information necessary to 
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find the required property of f and to store it into the output vector collection. 

[0085] In order to represent quantum gates, it is convenient to employ diagrams called quantum circuits. An exam- 
ple of quantum circuit is shown in Figure 1 1 . Each rectangle is associated with a matrix nxn, where n is the number of 
lines entering and leaving the rectangle. For example, the rectangle marked U F is associated to'matrix U F . Quantum 
5 circuits provide a high-level description of the gate. Using some transformation rules, one can compile them into the cor- 
responding gate-matrix. These rules are illustrated in Figures 12A-12F 

[0086] Jhe decoder block 806 interprets the basis vectors collected after the iterated execution of the quantum 
block 804. Decoding these vectors means to retranslate them into binary strings and interpreting them directly if they 
already contain the answer to the starting problem or using them, for instance, as coefficients vectors for some equation 
10 system, in order to get the searched solution. 

[0087] As an example of the above algorithm, it is instructive to develop a quantum gate representation of the Deut- 
sch-Jozsa's (DJ) algorithm. This gate is implemented according to the technique developed in connection with Figure 
8 to show the structure of matrix operator U F . 

[0088] The DJ algorithm is stated as: given a constant or balanced function £{0,1} n ->{0,1}, decide if f is constant 
15 or balanced. (This problem is very similar to Deutsch's problem, but it has been generalised to n>1.) 

A. First consider the encoder for the DJ algorithm for the special case where n=2. Thus: 



vxe{0,1}": f(x) = 1 



In this case the f map table is given by: F 



The encoder block takes f map table as input and encodes it into matrix operator U F , which acts inside of a 
complex Hilbert space. 

Function f is encoded into the injective function F, built according to the following statement: 

F:{0,1}" +1 ^{0,1}" +1 : F(x 0 ,x„y 0 )=(x 0 ,x v f(x 0 ,x,)<By 0 ) 

Then the F map table is: 



(*o.*i.yo) 


F(* 0 ,*i>yo) 


000 


001 


010 


011 ' 


100 


101 


110 


111 


001 


000 


011 


010 


101 


100 



12 



BNSDOCID: <EP 1083520A2_L> 



(continued) 



(*o>*i.yo) 


F{x b ,xi,y 0 ) 


111 


110 



F is encoded into U F map table using the rule: 

vte{0,1}" +1 : U F [T(f)]= x[F{t)) 
where t is the code map defined above. This means: 



1*0*1 yo> 




I000> 


1001 > 


101 0> 


1011 > 


I100> 


n oi > 


I110> 


im> 


1001 > 


I000> 


101 1> 


101 0> 


I101> 


I100> 


I111> 


I110> 



Starting from the map table of U F , calculate the corresponding matrix operator. This matrix is obtained using 
so the rule: 



u„ 


|000> 


|001> 


|010> 


|011> 


|100> 


|101> 


|110> 


llll> 


|000> 


0 


1 


0 


0 


0 


0 


0 


0 


|001> 


1 


0 


0 


0 


0 


0 


0 


0 


|010> 


0 


0 


0 


1 


0 


0 


0 


0 


|011> 


0 


0 


1 


0 


0 


0 


0 


0 


|100> 


0 


0 


0 


0 


0 




0 


0 


|101> 


0 


0 


0 


0 


1 


0 


0 


0 


|110> 


0 


0 


0 


0 


0 


0 


0 


1 


llll> 


0 


0 


0 


0 


0 


0 


1 


0 



So, U F is the following matrix: 

Using the matrix tensor product, U F can be written as: 

U F =I®I®C 

where ® is the tensor product, / is the identity matrix of order 2 and C is the NOT-matrix so defined: 
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HJ :] 



The matrix C flips a basis vector. In fact it transforms vector I0> into I1> and II > into I0>. If the matrix U F is 
applied to the tensor product of-three vectors of dimension 2, the resulting vector is the tensor product of the three 
vectors obtained applying matrix / to the first two input vectors and matrix C to the third. The structure of U F is such 
that the first two vectors in the input tensor product are preserved, whereas the third is flipped (action of C). This 
action corresponds to the constraints stated by U F map table. 
B. Now consider the case: 

n = 2 
vxe{0,1}": f(x) = 0 



In this case f map table is given by: 





f(x) 


00 


0 


01 


0 


10 


0 


11 


0 



The F map table is: 



(*o.*i.yo) 


F{x 0 ,x v y 0 ) 


000 


000 


010 


010 


100 


100 


110 


110 


001 


001 


011 


011 


101 


101 


111 


111 



The F map table is encoded into the Up map table: 
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t/plx 0 Xiyo> 






101 0> 


— ioioZ 


I100> 


I100> 


I110> 


I110> 


1001 > 


1001 > 


101 1> 


101 1> 


I101> 


I101> 


I111> 


I111> 



It is easy to transform this map table into a matrix. Every vector is preserved. Therefore the corresponding 
matrix is the identity matrix of order 2 3 . 





|000> 


|001> 


|010> 


|011> 


|100> 


|101> 


|110> 


llll> 


|000> 


1 


0 


0 


0 


0 


0 


0 


0 


|001> 


0 


1 


0 


0 


0 


0 


0 


0 


|010> 


0 


0 


1 


0 


0 


0 


0 


0 


|011> 


0 


0 


0 




0 


0 


0 


0 


|100> 


0 


0 


0 


0 


1 


0 


0 


0 


|101> 


0 


0 


0 


0 


0 


1 


0 


0 


|110> 


0 


0 


0 


0 


0 


0 


1 


0 


llll> 


0 


0 


0 


0 


0 


0 


0 


1 



Using matrix tensor product, this matrix can be written as: 
U F = I ® / ® / 

The structure of U F is such that all basis vectors of dimension 2 in the input tensor product evolve independ- 
ently. No vector controls any other vector. 

For example, consider the balanced function: 

n =2 

v(x 1 ,..,x n )e{0,1} n : f(x v ..,x n ) = x A ®..®x n 
In this case f map table is the following: 



X 


m 


00 


0 


01 


1 


10 


1 
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(continued) 



X 


fix) 


11 


0 



The following map table, calculated as described above, represents' the injective function F (where f is 
encoded into): 



w 





F(x 0 ,xi,y 0 ) 




(x 0 ,x v y 0 ) 


F(x 0 ,x u y 0 ) 


000 


000 




001 


001 


010 


011 




011 


010 


100 


101 • 




101 


100 


110 


110 




111 


111 



20 

Now encode F into U F map table: 



lx 0 x 1 y 0 > 


U F \XoXty 0 > 


I000> 


I000> 


101 0> 


101 1> 


I100> 


I101> 


I110> 


I110> 


1001 > 


I001> 


101 1> 


101 0> 


I101> 


I100> 


I111> 


I111> 



The matrix corresponding to U F is: 





|000> 


|001> 


j010> 


|011> 


\\00> 


|101> 


|110> 


llll> 


|000> 


1 


0 


0 


0 


0 


0 


0 


0 


|001> 


0 


1 


0 


0 


0 


0 


0 


0 


|010> 


0 


0 


0 


I 


0 


0 


0 


0 


|011> 


0 


0 


1 


0 


0 


0 


0 


0 


|100> 


0 


0 


0 


0 


0 


I 


0 


0 


|101> 


0 


0 


0 


0 


1 


0 


0 


0 


|110> 


0 


0 


0 


0 


0 


0 


1 


0 


llll> 


0 


0 


0 


0 


0 


0 


0 


1 
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This matrix cannot be written as the tensor product of smaller matrices. In fact, if written write it as a block 
matrix, it becomes: 
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|00> 


|01> 


|10> 


lll> 


|00> 


I 


0 


0 


0 


|01> 


0 


c 


0 


0 


|10> 


0 


0 


c 


0 


lll> 


0 


0 


0 


I 



This means that the matrix operator acting on the third vector in the input tensor product depends on the values 
of the first two vectors. If these vectors are I0> and I0>, for instance, the operator acting on the third vector is the 
identity matrix, if the first two vectors are I0> and 11 > then the evolution of the third is determined by matrix C. So, 
this operator creates entanglement, namely correlation among the vectors in the tensor product. 
C. Consider now a general function with n= 2. In this general case, the f map table is the following: 





W 


00 


f 00 


01 


'01 


10 


'10 


11 


'11 



with f}e{0,1}, /=00,01 ,10,11. If f is constant, then 3ye{0,1}vxe{0,1} 2 : f{x)=y. If f is balanced then 
30 \{ff. fj = 0}l=l{ ff. f-, = 1 }l The injective function F (where f is encoded) is represented by the following map table: 



1*0*1 yo> 


*>Fi*o*iyo> 


000 


OOfoo 


010 


01f 01 


100 


10f 10 


110 


11 '11 


001 


00 -n f 00 


011 


01 - f 0 i 


101 


io -./ 10 


111 





so Now encode F into the U F map table: 





u F i*o*iyo> 


I000> 


ioor 00 > 


101 0> 


I01f 0 i> 
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(continued) 





U F \x Q x^y 0 > 


I100> 


I10f 10 > 


I110> 


I11'11> 


1001 > 


I00 -. f 00 > 


101 1> 


101 -. f 0 T> 


I101> 


no-. f 10 > 


I111> 


in -. f u > 



The matrix corresponding to U F can be written as a block matrix with the following general form: 



V F 


|00> 


|01> 


|10> 


lll> 


|00> 


M» 


0 


0 


0 


|01> 


0 


Ki 


0 


0 


|10> 


0 


0 


M t o 


0 


lll> 


0 


0 


0 





where /W,=/ if f,=0 and M,=C if f}=1 , /'=00,01 ,1 0,1 1 . The structure of this matrix is such that, when the first 
two vectors are mapped into some other vectors, the null operator is applied to the third vector, generating a null 
probability amplitude for this transition. This means that the first two vectors are always left unchanged. On the con- 
trary, operators M,e{/, C] and they are applied to the third vector when the first two are mapped into themselves. If 
all M, coincide, operator U F encodes a constant function. Otherwise it encodes a non-constant function. If 
\{M-. Mj = /}I=I{W,-: Af; = C}\ then f is balanced. 
D. F or the general case n>0, the input function f map table is the following: 



X € {0,1}" 




0..0 


f o..o 


0..1 


f 0 ..1 






1..1 





with f)e {0,1}, ie {0,1}". If f is constant then 3ye {0,1}vxe{0,1}" : f(x)=y. If f is balanced then 
45 \\fj-. fj = 0}l=l{f,-: f/ = 1}l . The map table of the corresponding injective function F is: 



x e {0,1}"* 1 


F[x) 


0..0 


0-0 f 0 ..o 






1..10 


1-1 h.A 


0..01 


O-0-f-o..o 






1..11 


1-1 --'1..1 
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Now encode F into U F the map table: 



lx> 


U F \x> 


I0..00> 


I0..0 f 0 _ 0 > 






I1..10> 


11-1 f1..1> 


I0..01> 


I0..0 - f 0 0 > 






1..11> 


11. .1 ->f-l. A > 


0..1 





[0090] The matrix corresponding to U F can be written as a block matrix with the following general form: 





[0..0> 


|0..1> 




ll..l> 


\0..0> 


K..o 


0 


0 


0 


|0..1> 


0 


Mo.., 


0 


0 












ll..l> 


0 


0 


0 


A/,, 



where M,=/ if fj=0 and /W,=C if /y=1, ;'e{0,1} n . 
[0091] This matrix leaves the first n vectors unchanged and applies operator /W,e{/, C} to the last vector. If all M-, 
coincide with / or C, the matrix encodes a constant function and it can be written as "/ ®/ or "/®C. In this case no 
entanglement is generated. Otherwise, if the condition M - t = /}I=I{/W ( -: M,= C}\ is fulfilled, then f\s balanced and 
the operator creates correlation among vectors. 

[0092] The matrix U F , the output of the encoder, is now embedded into the quantum gate of the DJ algorithm. This 
gate is described using a quantum circuit shown in Figure 13. Using the identity rule of Figure 12C, the previous circuit 
compiles into the circuit shown in Figure 14. 

[0093] If f is constant and its value is 1 , the matrix operator U F can be written as "/®C. This means (as shown in 
Figure 12A) that U F can be decomposed into n+1 smaller operators acting concurrently on the n+1 vectors of dimen- 
sion 2 in the input tensor product. The resulting circuit representation is shown Figure 15. Using Figure 12B the sub- 
gate acting on every vector of dimension 2 in input, as shown in Figure 1 6. Every vector in input evolves independently 
from other vectors. This is because the operator U F does not create any correlation. So, the evolution of every input vec- 
tor can be analysed separately. This circuit can be written in a simpler way, observing that M' l=M as illustrated in Fig- 
ure 1 7. It can be shown that H 2 =/ , thus, the circuit can be redrawn as shown in Figure 1 8. 
[0094] Consider now the effect of the operators acting on every vector: 



. I0)-I1) 



so [0095] Using these results in connection with the operations shown in Figures 12D and 12C, gives the circuit rep- 
resentation shown in Figure 1 9. Thus, if f is constant with value 1 , the first n vectors are preserved. 
[0096] A similar analysis can be repeated for a constant function with value 0. In this situation U F can be written as 
n l®l and the final circuit is shown in Figure 20. In this case too, the first n input vectors are preserved. So, their output 
values after the quantum gate has acted are still I0>. 

55 [0097] The gate implementing the DJ algorithm in the general case is obtained operating on the circuit of Figure 14 
to yield thecircuit shown in Figure 21 , which is then evolved to the final circuit shown in Figure 22. 
[0098] If n=2, U F has the following form: 
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u. 


|00> 


|01> 


|10> 


lll> 


|00> 


M m 


0 


0 


0 


|01> 


0 


M m 


0 


0 


|10> 


0 


0 


Mo 


0 


lll> 


0 


0 


0 





w where M,e {/, C], (=00,01 ,10,1 1 . 

[0099] Given the quantum gate G=( 2 H®I) • U F * ( 2+1 H) , then: 



2 H<8/ 


I00> 


101 > 


I10> 


I11> 


I00> 


1/2 


1/2 


1/2 


1/2 


101 > 


1/2 


-1/2 


1/2 


-1/2 


I10> 


1/2 


1/2 


-1/2 


-1/2 


I11> 


1/2 


-1/2 


-1/2 


1/2 



U F - 3 H 


I00> 


101 > 


11 0> 


I11> 


I00> 


MoqH/2 


M 00 H/2 


MooH/2 


M 00 H/2 


101 > 


M 0A H/2 


-M 01 H/2 


M 01 H/2 


-M 01 H/2 


I10> 


M w H/2 


M 10 H/2 


-M 10 H/2 


-M 10 H/2 


I11> 


M U HI2 


-M U HI2 


-M„H/2 


M u H/2 



G 


I00> 


101 > 


11 0> 


I11> 


I00> 


(M 00 +/W 01 +/W 10 + 


(Moo-M 0 ^M w -M u )H/4 


(M 00 +M ov M w -M n )H/4 


fM 00 -/W 01 -M 10 +/W 11 ;H/4 


101 > 




^oo+/w 01 +yw 10 + 

M„)H/4. 


(Moo-M 0 yM 10 +M„)H/4 


^ 00 +/w 01 -/w 10 -/w 11 ;h/4 


I10> 


(M 00 +M or M AO -M„)H/4 


(Moo-M or M K +M^)H/4 


(M 00 -t./W 01 +/W 10 +/W 11 ;H/4 


(M m -M 0 ^M W -M^)HI4 
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(continued) 



G 


I00> 


101 > 


11 0> 


I11> 


I11> 


(M 00 -M 0VW +M„)H/4 


(M 00 +M ov Mi 0 -M n )H/4 


(M 00 -M w +M w -M n )H/4 


(M 00 +M 01 +/W 10 + Mn)H/4 


[0100] 


In the general case n>0, U F has the following form: 







|0..0> 


|0..1> 




ll..l> 


|0..0> 


K..o 


0 


0 


0 


|0..1> 


0 




0 


0 












ll..l> 


0 


0 


0 





where /W,e{/, C}, ie {0,1}". Given the quantum gate G=( n H<8>/) • U F '( n+1 H) , then: 





|0..0> 


!/> 


|1-.1> 


|0..0> 


H/2" a 






i> 


H/2" a 




. {-\y {lA) HI2 nn 




H/2* 1 


■(>•'> ' 


"(.!)(«...> 


11> 




JH/2 nn 


(i..i)^/2«« 



[0101 ] The above uses the binary string operator " • ", which represents the parity of the AND bit per bit between 
two strings. Given two binary strings x and y of length n, define: 

* • y = *i • yT ® x 2 • y 2 ©..© x„ • y„ 

where the symbol • used between two bits is interpreted as the logical AND operator. 
[0102] To prove that the matrix n+1 H really has the described form, it is useful to show that: 



45 [0103] The proof is by induction: 
, For n=1 : 



, , 1 (-1)<° 



[ 1h ]i,o = ^ = 
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„(/7-1)/2 -n/2 



1 nl 1 Ml (i1) ' (/,) 



[0104] The matrix n+1 H is obtained from n H by tensor product. Similarly, matrix n HS>l is calculated: 



"J9®/ 


|0..0> 




|i..i> 


|0..0> 


in* 1 




I/2* 2 


l> 




... (-xyjiw 1 


.. (-iy < i " , >//2« /j 


|ii> 


I/2" a 


... (Af- t) Jl/2» a . 


( . 1)( i..ihV..i, //2 „/2 





U F 


|0..0> 


|/> 


ll..l> 


30 


|0..0> 




M^I2 m 






\i> 


MiHll" a 


i-iy) 'Mfii2 m 


(-l)»' (, - ,) M/fl"/2" /2 


35 


Ii'.."i> 


M x . A HJ2 na 


... (-i) (IM) yM, 1 fl/2«' 3 





[0105] This is calculated only the first column of gate G since this operator is applied exclusively to input vector 
40 I0..01> and so only the first column is involved. 



G 


|0..0> 


|0..0> 


(Af 0 .. 0 +..+M,+..+M I .. l )/i r /2" 
(S j ,,,o,,«(-l) (, " ,> %W/2» ... 




\> 




|i..i> 



[0106] In the case of f constant, this means that all matrices M, are identical. This implies: 
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-l(Z y (-1)'' -/ M / )H = 0 



i since in this summation the number of +1 equals the number of -1 . Therefore, the input vector I0..01 > is mapped 

into a superposition of vectors I0..00> and I0..01>. If / is balanced, the number of M,.=/ equals the number of M,= C . 
This implies: 



— ri (ZjMj)HM) = 0 



[0108] This means that input vector I0..01>, in the case of balanced functions, can't be mapped by the quantum 
gate into a superposition containing vectors I0..00> or I0..01>. 

[0109] The quantum block terminates with measurement. It is possible to determine the possible outputs of meas- 
urement and their probabilities: 



Superposition of Basis Vectors 


Result of Measurement 


Before Measurement 


Vector 


Probability 


Constant functions: 


|0..00> 




CI0..01>=|0..0>O(oi|0>+a ( |l>) 


|0..01> 


IIor.il 2 


Balanced functions: 






G|0..01>=Z /)!(OtI1 „ . (0 00i0 „ 0 „ cti\i> 


Vje{0,l}«-{0»00, o..01}:|i> 


ll«,f 



[0110] The set A-B is given by all elements of A, unless those elements belonging to B too. This set is sometimes 
denoted as A/B. The quantum block is repeated only one time in Deutsch-Jozsa's algorithm. So, the final collection is 
made only by one vector. 

45 [0111] As in Deutsch's algorithm, when the final basis vector has been measured, it is desirable to interpret it in 
order to decide if f is constant or balanced. If the resulting vector is I0..0> the function is constant, otherwise it is bal- 
anced. In fact, the gate G produces a vector such that, when it is measured, only basis vectors I0..00> and I0..01 > have 
a non-null probability amplitude exclusively in the case f is constant. Besides, if f is balanced, these two vectors have 
null coefficients in the linear combination of basis vectors generated by G. In this way, the resulting vector is decoded 

so in order to answer Deutsch-Jozsa's problem: 



Resulting Vector after 
Measurement 


Answer 


I0..00> 


f is constant 
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(continued) 



Resulting Vector after 
Measurement 


Answer 


I0..01> 


f is constant 


otherwise 


f is balanced 



[0112] Graver's algorithm, is described here, is a variation on the Deutsch-Jozsa's algorithm. Graver's algorithm is 
10 stated as follows: 



Input 


given a function f: {0,1} n ->{0,1} such that 

axe {0,1}": (f(x)=lAVye {0,1}":x*y=>f(y)=0) 


Problem 


Find* 



[0113] In the Deutsch-Jozsa algorithm describes two classes of input functions and the problem is to decide what 
20 class the input function belonged to. In Graver's algorithm the problem is similar in form, even if it is harder because now 
one is dealing with 2" classes of input functions (each function of the kind described is considered as a class). In order 
to make the discussion more straightforward, the treatment that follows first considers a special function with n=2. Then 
the general case with n=2 is described, and finally the general case with n>0 is described. 
[0114] First consider the case where: 

25 

n=2 /(01)=1 
[0115] in this case the f map table is defined as: 

30 



X 


f(x) 


00 


0 


01 


1 


10 


0 


11 


0 



40 

[0116] The function f is encoded into injective function F, built as described in connection with Figure 8: 
F:{0,1} n+1 ^{0,1} n+1 : F(x 0 ,x v y 0 )=(x 0 ,x v f(x Q ,xJ@y 0 ) 
45 [0117] Then F map table is: 



(*o>*i.yo) 


F(*o.Jfi,yo) 


000 


000 


010 


011 


100 


100 


110 


110 


001 


001 


011 


010 
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(continued) 



(*o>*i>yo) 


F(x 0 ,x-„y 0 ) 


101 


101 


111 


111 



[01 18] Now encode F into the map table of U F : 
10 vse {0,1}" +1 :U f [t(s)]=t[F( S )] 

where i is the code map described in connection with Figure 8. 

15 



i*o*iyo> 


u F ix 0 *iyo> 


I000> 


I000> 


101 0> 


101 1> 


I100> 


I100> 


I110> 


I110> 


1001 > 


1001 > 


1011> 


101 1> 


I101> 


I101> 


I111> 


I111> 



[0119] From the map table of L/ F calculatethe corresponding matrix operator. This matrix is obtained using the rule: 
[U F ] ij =^^ U F \j) = \t) 

. 35 to give: 



[0120] The effect of this matrix is to leave unchanged the first and the second input basis vectors of the input tensor 
product, flipping the third one when the first vector is I0> and the second is 11 >. This agrees with the constraints on U F 
so stated above. 

[0121] Now take into consideration the more general case: 



n = 2 /(x) = 1 
55 [0122] The corresponding matrix operator is: 
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|00> 


|01> 


|10> 


|11> 


|00> 


Ho 


0 


0 


0 


|01> 


0 


M) t 


0 


0 


|10> 


0 


0 




0 


lll> 


0 


0 


0 


M n 



with M x =C^i*x:M^I. 

[0123] It is straightforward to generalize operator U F from the case n=2 to the case n>1 . The operator C on the 
main diagonal of the block matrix, is based on correspondence of the cell labeled by vector lx>, where x is the binary 
string having image one by Therefore: 





|00> 


|01> 


lll> 


|00> 


Moo 


0 


0 


|0l> 


0 


Ki 


0 












lll> 


0 


0 





25 with M x =CA\/i*x:Mi=t. 

[0124] The matrix U F , the output of the encoder, is embedded into the quantum gate. This gate is illustrated by the 
quantum circuit shown in Figure 23. The operator D n is called a diffusion matrix of order n and it is responsible for inter- 
ference in this algorithm. It plays the same role as QFT n in Shor's algorithm and of n H in Deutsch-Jozsa's and Simon's 
algorithms. This matrix is defined as: 

30 



D. 


|0..0> 


|0..1> 




jl..0> 


|1..1> 


|0..0> 


-1+1/2- 


1/2"-' 


1/2"- 1 




1/2"-' 


|0..l> 


1/2"-' 


-1+1/2-' 


1/2"-' 


1/2"" 1 


l/2»- 1 


\i> 


1/2"-' 


1/2"-' 


-1+1/2- 1 


1/2"-' 


1/2"" 1 


|1..0> 


1/2""' 


1/2"-' 


1/2"' 1 


-1+1/2"-' 


1/2"-' 


ll..l> 


1/2"-' 


1/2"-' 


1/2""' 


l/2 n -' 


-1+1/2-' 



45 [0125] Using the transformation shown in Figure 1 2C, the circuit of Figure 23 is compiled into the circuit of Figure 
24. 

[0126] For example, consider the case where U F has the following form: 
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[0127] Then calculate the quantum gate G=[(D 2 ®/) • U F ] h • ( 2+1 H) which in this case becomes: 



£ 


3 H 


I00> 


I01> 


11 0> 


I11> 




I00> 


H/2 


H/2 


H/2 


H/2 




101 > 


H/2 


-H/2 


H/2 


-H/2 


W 


I10> 


H/2 


H/2 


-H/2 


-H/2 




I11> 


H/2 


-H/2 


-H/2 


H/2 


15 














D 2 ®l 


I00> 


101 > 


11 0> 


I11> 




I00> 


-1/2 


1/2 


1/2 


1/2 


20 


I01> 


1/2 


-112 


1/2 


(a)/2 




I10> 


1/2 


1/2 


-1/2 


(b)/2 




I11> 


1/2 


1/2 


1/2 


-1/2 



U F - 3 H 


I00> 


I01> 


I10> 


I11> 


I00> 


H/2 


H/2 


H/2 


H/2 


101 > 


CH/2 


-CH/2 


CH/2 


-CH/2 


11 0> 


H/2 


H/2 


-H/2 


-H/2 


I11> 


H/2 


-H/2 


-H/2 


H/2 



[0128] Choosing /i=1 , yields: 



G 


I00> 


I01> 


11 0> 


I11> 


I00> 


(C+/)H/4 


(-C-/)H/4 


(C-3/)H/4 


(-C-/)H/4 


101 > 


(-C+3/)H/4 


(C+l)H/4 


(-C-/)H/4 


(C+l)H/4 


I10> 


(C+I)HI4 


(-C-0H/4 


(C+/)H/4 


(-C+3/)H/4 


I11> 


(C+l)H/4 


(-C+3/)H/4 


(C+/)H/4 


(-C-0H/4 



50 

[0129] In one example, the operator 3 H puts an initial canonical basis vector 1001 > into a superposition of all basis 
vectors with the same (real) coefficients in modulus, but with positive sign if the last vector is I0>, negative otherwise. 
The operator U F creates correlation: it flips the third vector if the first two vector are I0> and 11 >. Finally, D 2 ®/ produces 
interference: for every basis vector I xb*, it calculates its output probability amplitude 

55 
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by inverting its initial probability amplitude 
and summing the double'of the mean 

of the probability amplitude of all vectors in the form lx 0 *iJ4)>- In our example Oo=1/(4 • 2 1/2 ), a1= -1/(4 • 2 1/2 ). Take, for 
instance, basis vector I000>. Then 

a ' 0 oo=-«boo + 2«o=-l/(2-2 ,/2 )+2/(4-2" 2 )=0. 
[0130] In general, if n=2, U F has the following form: 



U F 


|00> 


|01> 


|10> 


lll> 


|00> 


M M 


0 


0 


0 


|01> 


0 


A*"o, 


0 


0 


|10> 


0 


0 


M l0 


0 


lll> 


0 


0 


0 


M n 



where M ; = CAVfex/Wp/fx/e :{0,1}") . 
[0131] The quantum gate G=(D 2 ® /) ♦ Up* ( 2+1 H) in this general case is: 



U F - 3 H 


I00> 


I01> 


I10> 


I11> 


I00> 


MqqH/2 


M 00 H/2 


MqoH/2 


M 00 H/2 


I01> 


M 01 tt/2 


-M m HI2 


/W 01 H/2 


-M 0 ,H/2 


I10> 


M 10 H/2 


M W HI2 


-M 10 H/2 


-M 10 H/2 


I11> 


M n H/2 


-M n Hl2 




M-hH/2 



G 


I00> 


I01> 


I10> 


I11> 


I00> 


(-/W 00 +/W 01 +M 10 + 


(-Moo-M 0 -\+M w - 


(-M 00 +/W 01 -/W 10 - 


(-M 00 -M or 




M^)H/4 


M U )HIA 


M^)HIA 


M w +M„)H/4 


101 > 




(Moo+Moi+Mn-M^H/A 




(M 00 +M or M w +M„)H/4 


11 0> 


(M 00 +M or M w +M u )H/4 


(/W 00 -M 01 -M 10 -M 11 )H/4 




(M 00 -M 01 +M 10 +/W 11 )H/4 


I11> 


(M 00 +M 0 - i +Mi 0 -M 1 ! )H/4 






{Moo-MovM-to-Mi^HM 



[0132] Now, consider the application of G to vector I001>: 
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(31001) = 1|00)®(-M 00 + M 01 + M 10 + M v )H\-\)+ 1|01>®(M 00 - M 01 + M 10 + M 11 )HI1)+ 
1|10>®(M 00 + M 01 -/W 1o + yW 11 )HI1>+l|11)®(/W O0 + M 01 + W 10 - M 11 )HI1) 

£ 

x=00: 
[0133] 

JO 

(31001)= 1|00)®(- C + 30HI1) + 1|01>®(C + 1|10)®(C + /)HI1> + 1|11)®(C + ) = lOO^i^iL^ 

x=01: 
15 [0134] 

GI001)= 1|00)®(C + /)HI1) -.- 1|01)®(-C + 3/)HM)+ 1|10>®(C + /)HI1) + 1|11)®(C + /)HI1> = I01>(^ll>j 
y=10: 

£0 

[0135] 

GI001)= 1|00>®(C + + 1|01)®(C + /)HI1>+ 1|10)®(-C + 3/)HI1) + 1|11)®(C + /)HI1> = hO^l^gUj 

25 x=11: 
[0136] 

GI001>= 1|00>®(C + /)HI1> + ±I01>®(C + /)HI1>+ 1|10>®(C + /)HI1) + 1|11)®(-C + 3/)HI1> = 11 D^^r) 

[0137] This means that if one measures the output vector and encodes it back the first two basis vectors of dimen- 
sion 2 in the resulting tensor product, the following results: 



X 


(a) esult 


(b) robability 


00 


00 


1 


01 


01 


1 


10 


10 


1 


11 


11 


1 



45 [0138] In the general case where n>0, U F has the following form: 



u. 


I0..O 


|0..1> 




!1..1> 


|0..0> 




0 


0 


0 


|0..1> 


0 




0 


0 












ll..l> 


0 


0 


0 





55 

where M x = CaV/#x:/W ( = l(xie :{0,1} n ) . . 
[0139] The quantum gate G=(D n ®l) h • U F '( n+1 H) is given by: 
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!0..o> 


[/> 


|1..1> 


|0..0> 


HI2 nn 


H/2" a 


HI2 na 


\i> 


77/2"* 


... (-iyJhw 2 


. (-iy l,A) H/2» a 


|ii> 


77/2"* 


... (-l) (, " ,) 777/2" /2 . 


'. (-i7'-'> 

•o» , >#'/2" /2 



Z>„®7 


|0..0> 


|0..1> 


|i> 


|1..0> 


|1..1> 


|0..0> 


-7+7/2-' 


7/2"- 1 


7/2- 


7/2- 


7/2- 


|0..1> 


7/2- 1 


-7+//2"- 1 


7/2"- 1 


7/2- 


7/2- 


\> 


//2 B -' 


//2 0 "' 


-7+7/2- 


7/2"-' 


7/2"-' 


|1..0> 


7/2 n -' 


7/2-' 


7/2-' 


-7+7/2"-* 


//2 n-. 


|1..1> 


112^ 


7/2"- 1 


7/2"-' 


7/2- 


-7+7/2- 





|0..0> 


!/> 


|1..1> 


|0..0> 






M 0 ..o77/2"' z 




Miff/2" 11 


{-\fMfll2 nn 


(-l) i(,J) M77/2" /J 


|i"i> 


M UA HI2 na 


... (-l) (, " ,)i A/ 1 .. I 77/2" /2 





[0140] Assuming, for example ri=1 , then: 



G 


|0..0> 


|0..0> 


(-M a . 0 +Iy. (0Jl « Mi/2»-')H/Z»" 


l«> 









[0141] Being M ; = C and \H*x:M-= I , this column may be written as: 
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G 


|0..0> 




|0..0> 

|i'..'i> 







10 

and so: 



G 


|0..0> 


|0..0> 

|i"..'i> 


{[-l+(2«-l V2"-']7+ C/l n ->}m*V4 
{(2 n -iy2»- l I+[-l+l/2 n ' l ]C}H/2" /2 
{[-l+(2"-l)/2"- l ]/+ C/2. n - x }H/2^ 2 



25 [0142] Now, consider application of the matrix operator {[-1+(2 n -1)/2 n " 1 ]/+ C/2"" 1 ) H/2 n/2 and matrix operator 
{ (2 n -i)/2 n - 1 / + [-i+i/ 2 n " 1 ]C}H/2 ra/2 to the vector 11 >, then: 

■*o [0143] This means: 

o M |a.oi>- 

- [(-.^)».-»)^-.^)o..0- + (m^) j )... + (-.^) 1 .,)]^ 



which can be written as a block vector: 

so 



G„ =1 I0..01> 


I0..0> 


[-1+(2 n -2)/2 n " 1 ]/2 n;2 HI1> 




[+1+(2"-2)/2"- 1 ]/2 n/2 WI1> 
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(continued) 


G A=1 I0..01> 


I1..1> 


[-1+(2"-2)/2''' 1 ]/2' ,/2 HI1> 



[0144] Now, application of the operator (D„®/) • U F to a vector in this form yields: 

10 





I0..0> 


aHI1> 




pHI1> 


I1..1> 


oHI1> 



where a and p are real number such that (2 "-1)a 2 +p 2 =1 . The result is: 
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[0145] This means that if one starts from vector G^IO-.OI >, which is in the form considered, and applies h times 
operator (D n ®l) • U F , the coefficients at time f are such that: 

(2"-1)a M -p M - 

a t =2 i-! — — - a M 

2 



(2"-1)a (1 -p M 
2" 



[0146] So, as (S increases, a decreases, 
is [0147] As in Deutsch-Jozsa's algorithm, when the output vectorfrom the Grover quantum gate has been measured, 
one must interpret it in order to find x. This step is relatively simple. In fact, it is sufficient to choose a large h in order to 
get the searched vector lx>IO> or lx>H> with probability near to 1 . After obtaining the vector, it is encoded back into 
binary values the first n basis vector in the resulting tensor product, obtaining the string x as final answer. 

20 Information Theory 

[0148] The evolution of a quantum algorithm from an information theory point of view is now disclosed. The com- 
plex vector entering the quantum gate is here considered as an information source both from the classical and the 
quantum level. The Shannon Entropy H Sh is a measure of classical information. Consider a complex vector of modulus 
25 1 in the Hilbert space 

Hii Qt ®...®«z a , 

so where 



35 has dimension 2 for every k, written as a complex linear combination of basis vectors: 



l V>= £ a / 1 / 2 .../ n l/ 1>® l/ 2>® - ® , 'n> 

/ 1 ,/ 2 ,.:.,;>{0,1} 



[0149] Then,' the Shannon Entropy of ly) with respect to the basis 



is defined as: 



H sh 0V» = - £ «" vV-// l0 9(H a / n / 2 .. / i 2 ) 

/ 1 ./ a ...../ l , 6 {0,1} 



55 where 
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M" 

is interpreted as the probability of measuring vector \i A i 2 ...i n ). 

[0150] The Von Neumann Entropy is used to measure the information stored in quantum correlation. Let p = l\|/>(yl 
be the density matrix associated to state and Tc {1 ,...,n}. Then define: 

Pt= ; 7r {1i .„ n) . T p 

where Tr {1 n} . T (...) is the partial trace operator. The Von Neumann Entropy of qubit j in ly> is defined as: 
S| V) (T)=-fr(p T logp T ) 
[0151] The following definitions are also useful: 

l^(h;k) = S, v) ({h}) + S lv) ({*}) - S iv) ({h,k}) 
l l>v) (h;k;l) = S lv) ({h)) + S^({k)) + S, v) ({/}) - S lv) ({h,k))-S^({h,l}) - S lv) ({k,l}) - S lv) ({h,k,l}) 
S {v) (h\k) = S^({h})-l i>0 (h;k) 
S l[V) (hlk;l) = S, v) ({/)}) - l iv) (h;k) - l^(h;k) + l^ } (h;k;l) 

for h,k,le 

[0152] These quantities can be represented in a Wenn diagram as shown in Figure 25. 

[0153] Measures of Entropy are different from most physical quantities. In quantum mechanics one has to distin- 
guish between observables and states. Observabies (like position, momentum, etc.) are mathematically described by 
self-adjoint operators in Hilbert space. States (which generally are mixed) are characterised by a density matrix p > 0, 
i.e. a Hermitian operator, with trace Tr(p) = 1 . The expectation value of an observable A in the state p is (A) = Tr(pA) . 
Entropy is not an observable. Thus, there does not exist an operator with the property that its expectation value in some 
state would be its entropy. It is rather a function of state. Due to the Jaynes relation between the information-theoretical 
and physical entropy via the Boltzmann's constant, k B , one can ascribe to any quantum object a certain value of its 
physical entropy 

S ph = k B H Sh- 

[0154] The classical limit S$ of the expression for the entropy can be justified mathematically on coherent states. 
The best thing one can do is to measure the probability of finding a particle in a state with minimum uncertainty centred 
around the classical values, i.e. in a coherent state. In the general case Sph > S, with the observation 




let \2) = W(z)\0) be a coherent state with expectation values of position, or momentum, q or p, respectively. In config- 
uration space, 10) is explicitly given by the wave function 




and W(z) is the unitary operator 
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with Q,P operators of position, or momentum, respectively. Now define the classical density distribution corresponding 
to the density matrix p by p(z):=(z\p\z) . 

[0155] For every function f{z) there exists at most one density matrix p that p (z) = f{z) and 
7rp = ^<zlplz>=1. 

The relation Sp h > S is true because, for s(x):= {- x In x(x > 0);0(x = 0)} due to concavity, S«zlplz)) > (zlplz) , 
is hence 

Spl = i~s(p(z)) >.\^(z\s(p)\z) = 7r(s(p)) = S(p). 

20 

More generally, for any convex (concave) function /, Tr{f(p))< (>) J^f f"(p(z)). By continuity of p(z), 



would imply S«zlplz» = <zlS(p)lz> forallz, i.e. regarding the strict concavity of S(«), every lz)mustbe an eigenvector 
of p, which is impossible . Hence, S$ > S. 

[0156] The classical entropy is not invariant under every unitary transformation, i.e. it cannot be said that 
S%(U' 9 U)=S%(p) 

for every U, but, rather, this relationship holds for a restricted class only. For instance if U = W(z 0 ) , then 

35 S°(W(-z 0 ) P W(z 0 ))=l^S«0\W(- Z )W(z 0 ) W(z 0 )W(z 0 )\0)) 

= ffs(p(z + z 0 )) 

= s>) 

40 

[0157] This argument also works for all unitary U such that UW(z) = W{z') times a phase factor provided that 
dz = dz' (canonical transformation). If \y) is a pure state (a unit vector), p = ly>(yl , then p(z) = l<ylz)l 2 , and 

4S < h (p) = -2/f Kvlz)l 2 lnl< V lz>l 2 

Inserting for ly) = lz 0 > , we obtain 

S p/l(P) = J — e =1- 

On the other hand, there exist pure states with arbitrary high classical entropy: it suffices to show that for every e > 0 
55 one can find unit vectors l\|i> such that <ylz> < e for all z. For them a well-known inequality telis us that Sp h (p) > - Ine 2 . 
We conjecture that the slate with minimal classical entropy are exactly given by the density matrices IzXzl, and conse- 
quently, Sp h (p) > 1 . In order that Sp h (p) be small, Suplp(z)l must be close to one, otherwise the inequality mentioned 
before gives a value too large for classical entropy. Now if Suplp(z)l is exactly equal to 1, then, by continuity, there is 
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some z 0 with p(z 0 ) = 1 , i.e. <z 0 lplz 0 ) = 1 . Since gpg < 1 this implies ||p|| = 1 and p(z 0 ) = lz 0 > ; on the other hand, 7r(p) = 
1 , hence p = lz 0 ){z 0 l because all other eigenvalues of p must be 0. 

[0158] The minimum-uncertainty coherent states for a harmonic-oscillator potential can be defined as those states 
that minimize the uncertainty relation of Heisenberg (leading to the equality in the uncertainty relations), subject to the 
added constraint that the ground state is a member of the set. They are considered to be as close as possible to the 
classical states. Beyond the harmonic-oscillator system, coherent states can be developed for quantum (Schrodinger) 
systems with general potentials and for general Lie symmetries. These states are called (general) minimum-uncertainty 
coherent states and (general) displacement-operator coherent states. There is also a different generalization of the 
coherent states of the harmonic-oscillator system. This is the concept of "squeezed" states. (Squeezing is a reduction 
of quadrature fluctuations below the level associated with a vacuum.) 

[0159] The even and odd coherent states for one-mode harmonic oscillator (Schrodinger cat states) are represent- 
atives of nofj-classical states. Schrodinger cat states have properties similar to those of the squeezed states, i.e. the 
squeezed vacuum state and the even coherent state contain Fock states with an even number of photons. 
[0160] In quantum mechanics, two non-commuting observables cannot be simultaneously measured with arbitrary 
precision. This fact, often called the Heisenberg uncertainty principle, is a fundamental restriction that is related neither 
to imperfection of the existing real-life measuring devices nor to the experimental errors of observation. It is rather the 
intrinsic property of the quantum states itself. The uncertainty principle provides (paradoxically enough) the only way to 
avoid many interpretation problems. The uncertainty principle specified for given pairs of observables finds its mathe- 
matical manifestation as the uncertainty relations. The first rigorous derivation of the uncertainty relation from the basic 
non-commuting observables (i.e., for the position and moment, 



It is a consequence of the properties of the Fourier transform that connects the wave functions of the system in the posi- 
tion and momentum representation. 

[0161] It is useful at this point to present quantum uncertainty relations (UR) in terms of entropy or information 
("entropic UR" -EUR). The usual "standard UR" (for standard deviations) 

(A^) 2 (A <p 6) 2 £ 1|<[4,B] j f I 2 ^-\ { {A,B) + ) V -2{A) V (B>/ 

(note that the second term in this inequality represents the covariance, or correlation, 
cov 9 (A,B):= l{q>l/4B + BA\q>) - «pl/Aiq»«piei<p> 

between the observables A and B in the state l<p)) is presented as an inequality of the entropic form 
S (A) + S (S) s s ^ Qr jn jnformatjon form y 4) + y B ) <■ y^ 6 ) as more adequate expressions for the uncer- 
tainty principle. Given two non-commuting observables, one can derive an uncertainty relation for them. The states that 
satisfy the equality sign in the inequality are called intelligent states. 

[0162] For example, consider any continuous parameter X and any Hermitian observable A(X) which is the gener- 
ator of the parametric evolution, then UR gives 
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is the parameter average of the observable uncertainty and 

10 K 

AX,= — (A..,-^) 
s o 

is the scaled displacement in the space of the conjugate variable of A. This generalized UR would hold for position- 
is momentum, phase-number or any combinations. For the case when initial and final states are orthogonal one knows 
that all states of the form 



are the only intelligent states which satisfy the equality 

However, these states do not satisfy the equality when the initial and final states are non-orthogonal. In this case, if 'the 
30 generator of the parametric evolution A can be split into two parts A 0 + 4 1 such that A 0 has a complex basis of normal- 
ised eigenvectors 

35 

which degenerate the spectrum fa 0 }, with / a set of quantum numbers and A, has matrix elements (A , ) „ = 0 = (A .,) » , 
and (A , ) /y - = (A 1 ) 9 = a 1 , then all states of the form: 

45 are intelligent states for non-orthogonal initial and final states. 

[0163] It is useful to compare various characterisations of "maximal information" and point out their connection with 
"minimum uncertainty." For clarity, the following treatment is restricted, without loss of generality, mainly to "simple" 
observables (defined on the smallest non-trivial Boolean algebra L = {0,a, ia,1}) . The quantity of interest is informa- 
tion with respect to the single effect E: l <f (E) = E <f ln(E v )+ E'^ ln(F ,,,), E'= I- E. Non-commutativity or incompati- 

so bility of (unsharp) properties E and F will, in general, exclude the possibility of measuring or preparing both of them 
simultaneously. In particular, if E = E Q (X),F = F p (Y) axe position and momentum spectral projections associated 
with bounded measurable sets X, Y, then E° (X)aE (Y) = 0 holds or, equivalently 

«plE Q (X)l<p>= 1 =» {<f\E p (Y)\<?)< 1 

55 

<(plE P (V-)lq»= 1 => «plE Q (X)l<p> < 1 
[0164] Thus "certain" position and momentum determinations exclude each other, and the question arises as to 
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what "degree of uncertainty" they can be "known" simultaneously. One can take any reasonable characterisation of 
maximal joint knowledge, or joint information. In this case the above-mentioned statement can be put into the following 
equivalent form: 



[0165] The "state of maximal information" can be defined through three values. The first expression £<p + F<p can be 
maximised and an explicit construction procedure for the corresponding "state of maximal information" is given below. 
The question of maxima for this quantity as well as for E^ • F 9 and for I^E) + I^F) for an arbitrary pair of effects, £ 
and F is studied here. In particular, it is shown that each quantity can be maximal only if there exist states which lead 
to minimal uncertainty product in UR. Furthermore, projections the maxima of /<p(E) + l v (F) (if they exist) coincide with 
those of one of the quantities £ 9 + F$ and £ 9 - F$ (£ v e {E,B),F^ e {F,F}). 

[0166] For maximal E 9 + F^ the variation of <<pl£l<p> + (<plFlq>) - X<cpl(p> must vanish which implies the following 
equations: (E + F)l<p) = (£ 9 + F 9 )lq» . Multiplying with £ or with F and taking the expectations yields: 

(A <p £) 2 = (A <p F) 2 = -(«pl£FI<p)-£ <p • F (p )=-cov (p (£,F), 

which leads to a minimal UR: (A^E) 2 • ( A 9 F) 2 = [cov 9 ( E,F)] 2 . 
[0167] Similarly, maximising the product • F 9 gives (F f £ + E.F)\<s>) = 2£ • F lq» and 

(A <p £) 2 F (p 2 = (A^F) 2 £ 9 2 = -E 9 • F^cov 9 (E,F) which leads again to a minimal UR, £^0*F f . 
[0168] Finally, maximal information sum l 9 (E) + yF) will be realised in states satisfying 
(InE - InE' ,,,)(£ - E' 9 )l<p> + (lnF 9 - lnF' 9 )(F - F 9 )lq» = 0 . Generally this equation contains all stationary points, e.g. 
the minimum 

E,-P f -F f -P f -J 

, or the joint eigenstates. Since one of the objects here is to find states of maximal information with respect to positive 
outcomes for E,F assume E v > \ and F v > \ . Then this equality implies: 



(<zE + F]<p) = (aF p +F v }<p),a = -XM * 0 

it) 



and 

45 

a(A 9 £) 2 = i(A ()) F) 2 = -cov(EF) 
which again gives rise to the minimal uncertainty product in UR. 

[0169] The three notions of maximal information are consistent in so far as they imply minimal uncertainty product. 
so For example, let E,F denote position and momentum spectral projections, respectively: £ = £° (X),F= F (Y) .The 
sum of probabilities E,, + F 9 has been shown to be maximal in the state <P = (p^ with 
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provided that X, Y are bounded measurable sets. Here a 2 is the maximal eigenvalue of the compact operator (FEF) 
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and gf 0 is the corresponding eigenvector satisfying 

FEF\g 0 ) = al\g 0 ),F\g 0 ) = lg 0 >,||lg 0 >|| 2 2 = 1. 

■5 It is clear from above description that cp min must be an eigenstate of (E + F). This can also be seen directly in the fol- 
lowing way. Introduce 

\f 0 )= a 0 1 £lgr 0 ),|y 2 2 = a-*{g o [FEF\g 0 ) = 1,E\f 0 ) = \f 0 ). 

io Then: 

EFE\f 0 ) = a 2 0 \f 0 ),\g 0 )=a^F\f 0 ) 
and (p min can be written in the symmetric form 

15 

l ^ = iEim lifo)+]9o)] 



so [0170] The value <p mln maximises all the three quantities (E 9 • F v ), (E^ + F^) and (y £) + l v (F)), and it minimises 
the uncertainty product A^E • A^F. Thus maximal information (minimal entropy) and minimal uncertainty can be 
achieved on intelligent coherent states and will again coincide. 

[0171] The following disclosure new describes the role of entropy changing in a quantum algorithm as information 
data flow processing and how the classical and quantum information amount changes in the dynamics of some quan- 
25 turn algorithms. The following qualitative axiomatic descriptions of dynamic evolution of information flow in a quantum 
algorithm are provided: 

(1) The information amount (information content) of a successful result increases while the quantum algorithm is in 
execution. 

30 (2) The quantity of information becomes the fitness function for recognition of successful results and introduces a 
measure of accuracy for successful results. In this case the principle of minimum of Classical/Quantum Entropy 
corresponds to recognition of successful results on intelligent output states of the quantum algorithm computation 
(3) If the classical entropy of the output vector is small, the degree of order for this output state is relatively larger, 
and the output of measurement process on intelligent states of a quantum algorithm gives the necessary informa- 

35 tion to solve the initial problem with success. 

[0172] These three information axioms mean that the algorithms can automatically guarantee convergence of infor- 
mation amount to a desired precision. This is used to provide robust and stable results for fault-tolerant computation. 
[0173] As an example of the use of entropy with quantum algorithms, consider the following examples from the 
40 Deutsch-Jozsa (DJ) algorithm with n = 3. Figures 27, 28, and 29 illustrate three input functions encoded into the follow- 
ing operators, respectively: 

[0174] The operator used in connection with Figure 27 is: 



7 


0 


0 


0 


0 


0 


0 


0" 


0 


c 


0 


0 


0 


0 


0 


0 


0 


0 


/ 


0 


0 


0 


0 


0 


0 


0 


0 


/ 


0 


0 


0 


0 


0 


0 


0 


0 


/ 


0 


0 


0 


0 


0 


0 


0 


0 


c 


0 


0 


0 


0 


0 


0 


0 


0 


c 


0 


,0 


0 


0 


0 


0 


0 


0 


c 



BNSDOCID: <EP 1083520A2J_> 



EP1 083 520 A2 



[0175] 



The operator used in connection with Rgure 28 is: 



5 




[0176] 



The operator used in connection with Rgure 29 is: 



w 




[0177] Figures 28-29 show the dynamics of the DJ algorithm, by monitoring the values for Shannon and Von Neu- 
mann Entropy at each step. Every step corresponds to the application of a quantum operator according to the general 
15 circuit diagram of Rgure 26. Examination of Figures 28-29 leads to the following conclusions about classical and quan- 
tum entropy changing after superposition, entanglement and interference have occurred. 

[0178] The input vector is a basis vector: the classical information of this state is 0. It is the tensor product of n basis 
vectors of dimension 2, so the Von Neumann entropy of every qubits composing it is also 0. 

[0179] The superposition operator 4 H increases the classical Shannon Entropy from its minimum value 0 to its 
20 maximum value 4, but does not change the entropy from the quantum Von Neumann Entropy point of view. 

[0180] The entanglement operator is a classical unitary operator, therefore it maps different basis vectors into dif- 
ferent basis vectors leaving the classical information of the system unchanged. However, it can create correlation 
among the different binary vectors in the tensor product describing the system state; this correlation is described by the 
Von Neumann Entropy of the different subparts of the system. The quantum information of the whole system is always 
25 0, even when the entanglement operator creates correlation, since the vector describing it is a pure state, whereas inner 
values for mutual information and conditional entropy can be positive or negative: they encode the quantum information 
necessary to decode the property being sought for the operator U F . The states of the system before and after the 
entanglement operator takes place cannot be distinguished from a classical information point of view, because the 
Shannon Entropy does not change. Only with a quantum information approach the difference between these two states' 
30 can be revealed. 

[0181] The interference operator leaves the quantum information picture unchanged, maintaining encoded the 
information used to identify U F as a constant or balanced operator. On the contrary, it decreases the classical entropy 
making the quantum information accessible; through the action of interference the vector acquires the minimum of clas- 
sical entropy: such a vector, according to the definition, is an intelligent state, because it represents a coherent output 

35 state of QA computation with minimum entropy uncertainty relation (EUR) as success result. 

[0182] Comparing Figures 27 and 28 shows that the entanglement operator in Figure 27 effectively creates quan- 
tum correlation among different parts of the system, whereas in Rgure 28 the general state is written as the tensor 
product of binary basis vectors and so no quantum correlation is involved. The interference operator in Figure 27 
reduces the classical entropy of 1 bit, whereas in Figure 28 it reduces the classical entropy of 3 bits. 

40 [0183] The presence of quantum correlation appears as the degree of resistance (immunity) of the system to 
change its classical entropy, and defines the internal degree of intelligent possibility of the quantum algorithm. 
[0184] The results in Figure 29 are similarto the results obtained in Figure 28. In Figure 29 the entanglement oper- 
ator creates no correlation. This is a common characteristic to all linear operators U F implementing a function f : {0,1 }" 
-» {0,1 } m such that f{x) = k • x or f(x) = ->{k - x) for some binary constant k, as shown in Figure 30. These func- 

45 tions, among the input set of balanced and constant functions, minimise to 0 the "gap" between the highest and lowest 
information values appearing in the Wenn Diagram of Figures 27-31 . 

[0185] The action of interference undergoes this property mapping of U F into an intelligent state revealing it. By 
contrast, other balanced functions are mapped into less intelligent states, that are higher classical entropy vectors. This 
means that it is a non-success result as it is shown in Rgure 31 . The Deutsch-Jozsa's algorithm undergoes the special 
so structure of its input set of functions from quantum information theory point of view. This structure is illustrated in Figure 
32. 

[0186] As a further example, Shor's algorithm is now considered from information theory viewpoint. The following 
operators implementing two input functions of period 2 and 4 respectively: 
[0187] The operator used in connection with Rgure 34 is (n = 3): 



40 
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L o c®c®c| 

[0188] The operator used in connection with Rgure 35 is (n - 2): 




"7®7 0 0 0 ' 

0 7®C 0 0 

0 0 C®7 0 

0 0 0 c®c 



[0189] Figures 34 and 35 illustrate the evolution of the quantum algorithm when applied with these operators. 
[0190] In Figure 34, the entanglement operator creates quantum correlation among vectors 3, 4 and 5: this corre- 
lation identifies the period of the input function. The interference operator preserves quantum correlation, but trans- 
poses it from vectors (3, 4, 5) to vectors (1, 4, 5). This transposing maintains the period of the input function encoded. 
But, it reduces the classical entropy from 3 bits to 2 bits, making it possible to access the period information generating 
an intelligent state, namely a state containing all the required quantum information but with minimum classical entropy 
as a qualitative measure of free energy. 

[0191] In Figure 35, the entanglement operator creates strong correlation among vectors 1 , 2, 3 and 4: this corre- 
lation identifies the input function that has maximal period (and so maximal entanglement). The interference operator 
preserves the correlation but does not decrease the classical entropy because entanglement is too great (the degree 
of resistance is too high). 

[0192] As shown in Figure 33, Shor's algorithm undergoes the special structure of its input space: periodical func- 
tions. Every function is characterised by its capacity to create quantum entanglement, which depends on its period. This 
structure is shown in Figure 33. 

[0193] Figures 36-38 illustrate information analysis for Graver's algorithm. The operator encoding the input function 
is: 

0 0 0" 
0 0 0 
0 0 0 
0 0 0 
0 0 0 
7 0 0 
0/0 
0 0/ 



"/oooo 

0 C 0 0 0 

0 0/00 

0 0 0 /0 

0 0 0 0 7 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 



[0194] Figure 36 shows a general iteration algorithm for information analysis of Graver's algorithm. Figures 37 and 
38 show two iterations of this algorithm. As shown in Figures 37 and 38, the entanglement operator in each iteration 
increases correlation among the different qubits. The interference operator in each iteration reduces the classical 
entropy but, as a side effect, it destroys part of the quantum correlation measured by the Von Neumann entropy. 
[0195] Graver's algorithm builds intelligent states in several iterations. Since every iteration first encodes the 
searched function by entanglement, but then partly destroys the encoded information by the interference operator, sev- 
eral iterations are needed in order to conceal both the need to have encoded information and the need to access it. 
[0196] The algorithms of Deutsch and Deutsch-Jozsa decision-making algorithms. The algorithms of Simon, Shor 
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and Grover are searching algorithms. The principle of minimum classical (quantum) entropy in output of a quantum 
algorithm means successful result are provided by intelligent output states. Decision-making algorthms provide recog- 
nition of intelligent output states with smaller or larger classical entropy than in the corresponding states of searching 
quantum algorithms. The searching quantum algorithms are based on a co-ordinated comparison of the minimum clas- 
5 steal entropy and the quantum entropy amount. The ability to co-ordinate these values characterises the intelligence of 
searching quantum algorithms. 

Claims 

10 i. A method for quantum algorithm design for quantum soft computing, comprising the steps of: encoding an input 
function 1 into a unitary matrix operator U F ; embedding the operator U F into a quantum gate G.said G being a uni- 
tary matrix; applying said quantum gate G to an initial vector to produce a basis vector; measuring said basis vec- 
tor; repeating said steps of applying and measuring k times, where 0 < k; and decoding said basis vectors, said 
. decoding including translating said basis vectors into an output vector. 

15 

2. The method of Claim 1, wherein said step of encoding comprises the steps of: transforming a map table of f into 
an injective function F; transforming said function F into a map table for said operator U F ; and transforming said 
map table for U F into said operator U F . 

20 3. The method of Claim 1 , further configured to minimize a Shannon entropy of said basis vectors. 

4. An intelligent control system comprising a quantum search algorithm configured to minimize Shannon entropy 
comprising: a genetic optimizer configured to construct one or more local solutions using a fitness function config- 
ured to minimize a rate of entropy production of a controlled plant; and a quantum search algorithm configured to 

25 search said local solutions to find a global solution using a fitness function configured to minimize Shannon entropy. 

5. The intelligent control system of Claim 4 wherein said global solution comprises weights for a fuzzy neural network. 

6. The intelligent control system of Claim 4 wherein said fuzzy neural network is configured to train a fuzzy controller, 
30 said fuzzy controller configured to provide control weights to a proportional-integral-differential controller, said pro- 
portional-integral-differential controller configured to control said controlled plant. 

7. A method for evolving a quantum search algorithm comprising the step of selecting a fitness function configured to 
minimize Shannon entropy. 

35 

8. A method for evolving a quantum search algorithm comprising the steps of: minimizing Heisenberg uncertainty; 
and minimizing Shannon entropy. 

9. A method for evolving a quantum algorithm comprising the steps of: applying an entanglement operator to create 
40 a plurality of correlated state vectors from a plurality of input state vectors; and applying an interference operator to 

said correlated state vectors to generate an intelligent state vector, said intelligent state vector having less classical 
entropy than said correlated state vectors. 

10. A method for global optimization to improve a quality of a sub-optimal solution comprising the steps of: applying a 
45 first transformation to an initial state to produce a coherent superposition of basis states; applying a second trans- 
formation to said coherent superposition using a reversible transformation to produce coherent output states; 
applying a third transformation to said coherent output states to produce an interference of output states; and 
selecting a global solution from said interference of output states. 

so 11. The method of Claim 10, wherein said first transformation is a Hadamard rotation. 

12. The method of Claim 10, wherein each of said basis states is represented using qubits. 

13. The method of Claim 10,. wherein said second transformation is a solution to Shrodinger's equation. 

55 

14. The method of Claim 10, wherein said third transformation is a quantum fast Fourier transform. 

15. The method of Claim 10, wherein said step of selecting comprises finding a maximum probability. 
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16. The method of Claim 10, wherein said superposition of input states comprises a collection of local solutions to a 
global fitness function. 

17. An apparatus for quantum soft computing, comprising: an encoder module configured to encode an input function 
s into a unitary matrix operator; an embed module configured to embed said unitary matrix operator into a quantum 

gate; a processing module configured to apply said quantum gate to an initial vector to produce a basis vector; 'a 
measurement module configured to measure said basis vector; and a decoder configured to decode said basis 
vector and translate said basis vector into an output vector. 

w 18. The apparatus of Claim 1 7, wherein said encoder comprises: a first transform module for transforming a map table 
of said input function into an injective function; a second transform module for transforming said injective function 
into a map table for said unitary matrix operator; and a third transform module for transforming said map table for 
said unitary matrix operator into said unitary matrix operator. 

15 19. The apparatus of Claim 17, further comprising a module to minimize a Shannon entropy of said basis vector. 

20. A method for intelligent control comprising a quantum search algorithm configured to minimize Shannon entropy 
comprising: optimizing one or more local solutions using a fitness function configured to minimize a rate of entropy 
production of a controlled plant; and searching, using a quantum search algorithm to search said local solutions to 

20 find a global solution using a fitness function to minimize Shannon entropy. 

21. The method of Claim 20 wherein said global solution comprises weights for a fuzzy neural network. 

22. The method of Claim 21 further comprising: training a fuzzy controller; providing control weights from said fuzzy 
25 controller to a proportional-integral-differential controller; and using said proportional-integral-differential controller 

to control said controlled plant. 

23. An apparatus for evolving a quantum search algorithm comprising: a module for calculating quantum entropy; a 
module for calculating classical entropy; and a module for searching a solution space for a solution that reduces 

30 both quantum entropy and classical entropy. 

24. An apparatus for evolving a quantum algorithm comprising: a first module for applying an entanglement operator to 
create a plurality of correlated state vectors from a plurality of input state vectors; and a second module for applying 
an interference operator to said correlated state vectors to generate an intelligent state vector, said intelligent state 

35 vector having less classical entropy than said correlated state vectors. 

25. A global optimizer to improve a quality of a sub-optimal solution, said optimizer comprising a computer software 
loaded into a memory, said software comprising: a first module for applying a first transformation to an initial state 
to produce a coherent superposition of basis states; a second module for applying a second transformation to said 

40 coherent superposition using a reversible transformation to produce one or more coherent output states; a third 
module for applying a third transformation to said one or more coherent output states to produce an interference of 
output states; and a fourth module for selecting a global solution from said interference of output states. 

26. The optimizer of Claim 25, wherein said first transformation is a Hadamard rotation. 

45 

27. The optimizer of Claim 25, wherein each of said basis states is represented using qubits. 

28. The optimizer of Claim 25, wherein said second transformation is based on a solution to Shrodinger's equation. 
so 29. The optimizer of Claim 25, wherein said third transformation is a quantum fast Fourier transform. 

30. The optimizer of Claim 25, wherein said fourth module is configured to find a maximum probability. 

31. The optimizer of Claim 25, wherein said superposition of input states comprises a collection of local solutions to a 
55 global fitness function. 
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An apparatus for quantum soft computing, comprising: means for encoding an input function f into a unitary matrix 
operator U F ; means for embedding the operator L/ F into a quantum gate G; means for applying said quantum gate 
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G to one or more initial vectors to produce one or more basis vector and measuring said basis vectors to produce 
measured vectors; and decoding said measured vectors into an output vector. 

33. The apparatus of Claim 32, wherein said means for encoding transforms a map table of f into an injective function 
5 F, transforms said function F into a map table for said operator U F and transforms said map table for U F into said 

operator U F 

34. The apparatus of Claim 32, further comprising means for minimizing an entropy of said basis vectors. 

10 35. An intelligent control system comprising a quantum search algorithm configured to minimize Shannon entropy 
comprising: means for optimizing one or more local solutions using a fitness function configured to minimize a rate 
of entropy production of a controlled plant; and means for quantum searching said local solutions to find a global 
solution using a fitness function configured to minimize Shannon entropy. 

75 36. The intell.jent control system of Claim 35 wherein said global solution comprises weights for a fuzzy neural net- 
work. 

37. The intelligent control system of Claim 35 further comprising means for training a fuzzy controller using said global 
solution, said fuzzy controller providing weights for a proportional-integral-differential controller. 

20 

38. An apparatus for evolving a quantum search algorithm comprising: a genetic optimizer based on a fitness function; 
and means for selecting said fitness function to minimize classical entropy and quantum entropy. 

39. An apparatus for evolving a quantum search algorithm comprising: means for minimizing a quantum entropy; and 
25 means for minimizing a classical entropy. 

40. An apparatus for evolving a quantum algorithm comprising: means for applying an entanglement operator to create 
a plurality of correlated state vectors from a plurality of input state vectors; and means for applying an interference 
operator to said correlated state vectors to generate an intelligent state vector, said intelligent state vector having 

30 less classical entropy than said correlated state vectors. 

41. An apparatus for global optimization to improve a quality of a suboptimal solution comprising: means for applying 
a first transformation to an initial state to produce a coherent superposition of basis states; means for applying a 
second transformation to said coherent superposition using a reversible transformation to produce coherent output 

35 states; means for applying a third transformation to said coherent output states to produce an interference of output 
states; and means for selecting a global solution from said interference of output states. 

42. The apparatus of Claim 41 , wherein said second transformation is a solution to Shrodinger's equation. 
40 43. The apparatus of Claim 41 , wherein said third transformation is a quantum fast Fourier transform. 

44. The apparatus of Claim 41 , wherein said means for selecting finds a maximum probability. 

45. The apparatus of Claim 41 , wherein said superposition of input states comprises a collection of local solutions to a 
45 global fitness function. 
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